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ABSTRACT 



The aim of this dissertation is to present a new method of 
engineering analysis and design for complex control systems , This 
method is the time domain infinite matrix method * The formulation of 
the infinite matrix follows from the convolution summation of sampled 
data systems. The mathematical basis of the time domain mtrix for- 
mulation is presented in a discussion of the applicable concepts of 
infinite matrices and sequence spaces. This method of analysis and 
design is applicable to both continuous data and sampled data systems. 
For continuous systems it is necessary to introduce a fictitious sam- 
pler and hold of sufficient sampling rate to effect an accurate ap- 
proximation , 

The time domain matrix method is presented and illustrated as a 
method of analysis and design of linear, nonlinear, and time varying 
systems of the continuous or sampled data class. Sampled data, time 
varying systems may not be conveniently investigated by any other 
existing method. Furthermore the investigation of nonlinear systems 
is greatly simplified by the time domin approach, Hiltiloop sys- 
tems may be treated with ease emd the signals at intermediate points 
throughout the loops are readily available. Also, systems with 
multiple nonlinearities may be investigated, for which there is not 
a presently available method of analysis and design. 

Two methods of design of a discrete compensator for a sampled 
data system are presented. These methods are accomplished directly 
in the time domain and allow for a compromise of specifications in 
the time domain. Also the response between sampling instants is 

ii 



accounted for in one of the two design procedures. 

The time domain matrix method may be readily programmed on a 
digital computer and therefore provides a rapid analysis and design 
technique . 



iii 



TABLE OF CONTENTS 



Chapter Page 

I INTRODUCTION 1 

II THE TIl^ DOMAIN MATRIX 

2.1 Introduction 5 

2.2 A Sampled Data System and the z Transform 5 

2.3 The Convolution Summation 7 

2.4 The Evaluation of the System Transfer 

^fe,t^ix of A Transfer Function 11 

2.5 Response Between the Sampling Instants 12 

2.6 The Formulation of the Ih,trix Equation 

for Closed Loop Systems 15 

2.7 The Solution for the Response of A 

Simple System 17 

III THE MATHEMTICAL THEORY OF INFINITE MATRICES 

3.1 The Mathematical Theory 20 

3.2 Convergence and Stability in the Sequence 

Space 22 



IV ANALYSIS OF CLOSED LOOP CONTROL SYSTE^6 

4.1 A Msthod of Evaluating the Response of 
A Closed Loop System Without Inversion 

of Mitrices 27 

4.2 Evaluation of the Response of A Closed 

Loop System With Other Than Unity 
Feedback 32 

4.3 Analjrsis of Jiiltiloop Control Systems 35 

4.4 Analysis of Time Varying Control Systems 39 

4.5 Analysis of Control Systems with Non- 

linear Components 47 

4.6 Analysis of Continuous Control Systems 53 

V INTRODUCTION TO THE DESIGN OF CLOSED LOOP 
CONTROL SYSTE.^B 

5 .1 Introduction 57 

5.2 The Design of the System Utilizing Time 

Varying Gain Compensation 59 

5.3 Application of the Time Domain ^fetrix 
tfethod to the Design of Various *^pes 

of Control Systems 62 



iv 



TABLE OF CONTENTS 



Page 



Chapter 



VI 


THE 


ANALYSIS AND DESIGN OF ADAPTIVE CONTROL 






SYSTEMS 






60 I 


Introduction 


72 




6.? 


The Use of A M^del in an Adaptive System 


73 




6o3 


An Adaptive Gain System 


78 




6.4 


An Adaptive Sampling Frequency System 


80 




6.5 


The Use of the Time Domain l^trix Msthod 
for Analysis and Design of Adaptive 
Systems 


85 


VII 


THE 


DESIGN OF DISCRI’TE DATA COMPENSATORS 






7.1 


Litroduction 


37 




7.2 


The Physical Realizability of Discrete 
Data Compensators 


88 




7.3 


The Sensitivity Matrix of the Sampled 
Data System with A Discrete Compensator 


90 




7.4 


Design of A Closed Loop Discrete Com- 
pensator 


91 




7.5 


Design of the Closed Loop Compensator 
by Time Domain Evaluation 


96 




7,6 


The General Characteristics of the Closed 
Loop Discrete Compensator 


99 




7.7 


Time Domain Design of A Discrete Com- 
pensator 


102 


VIII 


CONCLUSIONS 






8.1 


Summary of Results 


125 




8.2 


Further Conclusions 


126 




8,3 


Future Work 


128 




REFERENCES 


130 




BIBLOGRAPHY 


131 




APPENDICES 




A 


THE 


SYSTEM TRANSFER MATRIX 






A.l 


Determination of the System f&trix Usirig 





the z transform 132 



V 



TABLE OF CONTENTS 



Chapter 


A.2 


The Use of the Impulse Response to 


Page 






Evaluate the System Transfer I-feitrix 


134 




A,3 


Evaluation of the System Ti'ansfer >ktrix 
for Tnte^sample Response Using the 
Mfcdified z Transform 


135 




Ao4 


Evaluation of the System Transfer ^ktrix 
for Intersample Response Using the 
Lupulse Response 


137 




A. 5 


The System fktrix for A Type I Second 
Order System With A Hold 


138 




A. 6 


The System ^ktrix for A Type II Second 
Order System With A Hold Network 


139 


B 


A METHOD OF INVERTING A LOWER TRIANGULAR 
MATRIX 


U4 



Vi 



LIST OF ILLUSTRATIONS 



Figure 




Page 


1-1 


Sampled Data Control System Block Diagram 


2 


2-1 


Open Loop Sampled System 


6 


2-2 


An Open Loop Continuous System 


6 


2-3 


Open Loop Sampled System With Hold 


13 


2-4 


The Ihtersample Response 


13 


2-5 


First Order Sampled System 


16 


2-6 


Two Block Open Loop System 


16 


2-7 


Closed Loop Sampled System 


19 


2-8 


Closed Loop System With A Step Input 


19 


3-1 


Closed Loop Sampled Data System 


23 


4-1 


Error Sampled System 


27 


4-2 


Two Sampler Feedback System 


33 


4-3 


Error Sampled System 


33 


4-4 


Three Sampler Single Loop System 


36 


4-5 


Miltiloop Sampled System 


36 


4—6 


A Hiltiloop Sampled System 


38 


4-7 


Sampled Open Loop System 


40 


4—8 


Time Varying Closed Loop System 


41 


4-9 


System With Two Time Varying Components 


42 


4-10 


Time Varying System With A Feedback Component 


42 


4-11 


Open Loop Time Varjring System 


43 


4-12 


A Time Varying Gain 


44 


4-13 


System With A Time Varying Gain 


45 


4-14 


Time Varying Gain Example 


46 


4-15 


Single Loop Nonlinear System 


47 


4-16 


Saturating Amplifier Characteristic 


47 


4-17 


Saturating Amplifier Characteristic Example 


49 


4-18 


Hysteresis Characteristic 


50 


4-19 


System With Two Nonlinear Components 


51 


4-20 


Continuous Feedback Control Astern 


54 


4-21 


Location of The Fictitious Sampler 


55 


4-22 


The Approximation of A Time Function 


55 


5-1 


Unity Feedback Sampled Data Control System 


59 


5-2 


The Compensated Control System 


60 


5-3 


Step Response of Tfecoinpensated and Compensated 






System 


63 


5-4 


Continuous Data Compensation 


64 


5-5 


The Digital Compensator in A Sampled Data. 






System 


65 


5-6 


Control System With A Disturbance Input 


67 


5-7 


Continuous Control System With Fictitious 






Sampling 


70 


6-1 


The Use of A Msdel to Shape the Input Signal 


74 


6-2 


The Use of A Madel in A Conditional Feedback 






System 


74 



vii 



LIST OP ILLUSTRATIONS 



Figure Page 

6-3 A Sampled Data Conditional Feedback System 75 

6-4 Output Response for A Conditional Feedback 

System 77 

6-5 Adaptive System With A Variable Gain 78 

6-6 Compensated and Uhcoispensated System Response 81 

6- 7 Variable Sampling Rate System 82 

7- 1 A Sampled Data System With A Discrete Data 

Compensator 88 

7-2 The Sensitivity of the Compensated and Uncom- 
pensated System 92 

7-3 Error Sampled System With A Discrete Compen- 
sator 93 

7-4 An Open Loop Discrete Compensator 93 

7-5 Response of the System to A Step Input 97 

7-6 Response of the System to A Ramp Input 98 

7-7 System Response to A Unit Step Input 100 

7-8 The Stable Sampled Dapulse Response of the 

Discrete Compensator 103 

7-9 CoE^ensated Sampled Data Control System 105 

7-10 Compensated and Unconpensated System Response 109 

7-11 The Relation of the Intersample Response and 

the Sampled Response 115 

7-12 Compensated and Uncompensated System Response 116 

7-13 Compensated and Reduced Gain System Response 119 

7-14 The Relation of Intersample Response and the 

Sampled Response 120 

7-15 System Response to A Ramp Baput 122 

A-1 Open Loop Sampled System 132 

A-2 Itelayed Sampling 135 

A-3 The Impulse Response for Three Systems 143 



viii 



LIST OF TABIES 



Table 




Page 


5«1 


Rise Time arid Overshoot of the Coispensated 
and Uncompensated System 


62 


6„1 


Response with Fixed and Variable Sampling 
Period 




7,1 


Indices of Response for the Compensated arid 
Uncompensated Systems 


110 


A,1 


Element Values for the System Transfer 
l^trix for Various Systems 


140 



ix 



TABLE OF SYMBOLS AND ABBREVIATIONS 



Symbol 


Description 


s 


llie Laplace complex variable 


z 


sT 

z-transform variable (-e ) 


T 


The sampling period or interval 


r(t), R(s), R(z) 


Input Signal 


C(t), C(s), C(z) 


Output Signal 


D(z) 


Transfer function of discrete compensator 


G(s), G(z) 


Transfer function of the plant or system 


k,n 


An interger 


K 


Gain 


[G] 


System Transfer l^trix 


W 

r] 


Discrete Compensator ^ktrix 
Biput sequence column matrix 


GJ 

r (t) , R* (s) 


Output sequence column matrix 
Sampled Input Signal 

A sampler 



X 



CHAPTER I 
INTRODUCTION 

During the last two decades, feedback control systems have beocme 
increasingly important to our technological civilization, and have par- 
ticular] y contributed to national defense. Simultaneously, there has 
been increased interest and effort in the investigation of automatic 
control systems id.th respect to their analysis and design with various 
excitation signals , 

Hie I.R.E. defines a feedback control system as a control system 
comprising one or more feedback control loops, which combines functions 
of the controlled signals with functions of the commands to tend to 
maintain prescribed relationships between the commands and the con- 
trolled signals. Feedback control systems are a large class of sys- 
tems which include many subclasses of which a few are linear systems, 
nonlinear systems, multivariable systems, tin^ varying systems, sampled 
data systems , and adaptive systems , The active and dynamic elements 
are not limited and may be electronic, electromechanical, hydraulic, 
or pneumatic. 

An important subclass of control systems, sampled-data control 
systems, is a dynamical system which operates \d.th sampled or quantized 
information. That is, the information is present as a sequence of dis- 
crete numbers in time, in contrast to continuous data systems for which 
the controlling iiiformation is monitored continuously in time , Ordi- 
narily, the information is carried in the amplitude of the samples and 
may be considered pulse amplitude modulation. The block diagram of a 
typical sampled-data control system is shown in figure 1-1, where the 
sampled-controller may be a special or general purpose digital computer 
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Figure 1-1 „ Sampled Data Control System Block^Diagraffi 
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Kie wide use of digital computers today has further stimulated inves- 
tigation of sampled-data systems. The seemingly unlimited possibilities 
of the use of a digital computer as an active element allow the possi- 
bility of considering these systems as a beginning towards the realiza- 
tion of the much discussed adaptive and learning control systems which 
may imitate the human brain and nervous system. 

There exist at present three major methods for the analysis and 
design of sampled-data systems as follows ° 

1, Difference equation approach 

2, Frequency response methods 

3, Z transform and modified transforms and the use of the root 
locus , 

All these methods assume that sampling occurs instantaneously or there- 
fore, that the pulses are of negligible width. The first method uses 
classical methods of difference equations and yields a solution only at 
the sampling instants , Frequency response methods are an attempt to 
extend concepts from continuous systems and suffer from the attendant 
disadvantages , The z transform uses a complex-variable transformation 
and determines the performance of the system by location of roots on 
loci in the complex z plane, and by inverse transformation. 

It is the purpose of this dissertation to present a method of en- 
gineering analysis and design for complex control systems. This method 
is the use of infinite matrices in the time domain. Analysis and de- 
sign for this method takes place directly in the time domain and avoids 
the necessity of transformations in complex-variables. This method has 
proven to be potentially useful for the investigation of a variety of 
systems , 
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Time docain analysis and design has the important advantage of 
affording investigation directly in the domain of interest and direct 
evaluation of performance , therefore avoiding use of correlation 
theorems which are complex and may be inaccurate „ This method avoids 
the difficulty of solution that is present for higher order systems 
using the z-transformationc In fact, the amount of work necessary for 
investigation is approximately the same for a first order as for any 
n^^ order system. 

This method can be applied to many classes of systems which either 
do not lend themselves or are not possibly investigated by frequency or 
z-transform techniques. Therefore, beyond linear system investigation, 
one may analyze and design nonlinear and time-varying systems. Also 
the investigation of continuous systems is made possible through the 
introduction of a fictitious sampler and hold of suitably high sampling 
rate. The error introduced by this approximation of a continuous sys- 
tem by a sampled-data system can be reduced to a negligible amount with 
an attendant increase in labor of calculation. 

This dissertation is presented in three parts. The first part 
comprising Chapters 2 and 3 presents the mathematical background and 
formulation of the time domain mtrix equations. The second part com- 
prised of Chapters 4 and 5 presents the analysis and design methods and 
verification of the theory by application to various t 3 rpes of systems . 
The third part comprised of Chapters 6, 7, and 8 is concerned with the 
design and realizibility of the digital compensator in the time domain, 
and the final conclusions and possibilities for future investigation. 
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CHAPTER 2 



THE TIME DOMAIN MATRIX 

2-1 Introduction 

The method of investigation of sampled-data control systems most 
commonly used today is the z-transf ormation , as is correspondingly the 
s-transform for continuous systems . The z-transformation converts the 
difference equations to a set of simultaneous algebraic equations which 
may be solved for the unknown variable and then by means of the inverse 
transform yields the response in the time domain. Formulation of the 
problem directly in the time domain allows one to avoid this transfor- 
mation and inverse transformation, 

2-2 A Sampled-data System and the z Transform 

A block diagram of a simple open loop sampled-data system is shown 
in figure 2-1, The definition of the z-transform of x (t) is 

oo 

X(z) = X^ (~— In z) = x(nT)z“^ (2-l) 

n=0 

-sT 

where z=e =u + jv 

and T = sampling period 

Then it can be shown that the output at the sampling instants Y(z) is^ 

Y(z) = G(z) X(z) (2-2) 

where G(z) is the z transform of G(s), 

In order to obtain the response between the sampling instants 

2 

Jury introduced the modified z-transform where the response between 
the samples is^ 

Y(zjm) = G(z,ra) X(z) (2-3) 
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y(t) 



Y(s) 



l._ _ Y*(s) 

Fictitious Sampler 

Figure 2-1 » Open Loop Sampled System 



X(s) 

x(t) 



Syncronized Samplers 
with period T 




x(t) 



G(s) 



Y(t) 



Figure 2-2 <, Open Loop Continuous System 



In order to obtain the response in the time domain, the inverse 



transformation must be utilized of which one form is ° 




G(z)z*^ ^ dz 



(2-A) 



where is the path of integration in the z-plane 
which encloses all the stoigularities of the integrand,, 



It is found that the usefulness of the transform lies with the use of 
the root loci on the z-plane „ However, one desires to obtain the res- 
ponse directly in the time domain and that will be the subject of the 
next section, 

2-3 The Convolution Summation 

The time response of a continuous system G(s) as shown in figure 
2-2 to a continuous input x(t) is given by the convolution integral” 

t 

y(t) = \ g(t-X) x(X) dX (2“5 

-00 

where X = duran^r variable time delay 



This method of obtaining the response is usually avoided in favor of the 
Laplace transform due to the difficulty of evaluation of the integral. 
However, for the sampled data system as shown in figure 2-1, the 
input signal may be written^ 



and g(t“X) = the delayed impulse response of the system G(s), 



OO 




x(nT) &(t-nT) = 



(2-6) 



n=0 



c 



where |^x] is a column 



matrix of n elements 
where n extends to infinity 



X 



n 
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Now, in contrast to equation 2-5 » one may write an equation for the con- 
tinuous output response as a convolution summation" 

n 



y(t) = 




k=0 



g(nT -kT) x(kT) 



(2-7) 



n 



n 



or 



- ° ^ gkVk 

k=0 k=0 

where is often referred to as the weighting sequence. 

The values at the sampling instants of the impulse response, g^, are 
related to g(t) the impulse response as^ 



g^ = g(nT) = g(t) 



(2-9) 



t = nT 

Therefore, vdien the input is an ideal impulse of unity height, theni 
'"l for k = n 



^n-k ■ ^ 



k / n 



and y^ = g^ the weighting sequence, 



If equation 2-8 is expanded, one obtains! 



n 


= 0 


7o 


- go-0 


Xo 






n 


= 1 




= gl-0 


Xo 


+ g1-1 


Xi 


n 


= 2 


72 


” S2-0 


Xo 


+ 82-1 


Xl 



S2-2 ^2 



(2-10) 



^ ^ Vji ~ Sn—o ^0 Sn— 11 ^ 8 n— 2 

* ^n-n 

And equation 2-10 may be vrritten as ° 



O o O 



n 



8 



Vo - So 



Xo 



7i 


= 


Xq 


+ go X, 


(2-11) 


7z 

0 


= Sz 
0 


Xo 


+ gl Xi + go X 2 




0 

o 

^n 


o 

« 

= ^n 


Xq 


+ gn-1 X, + . . . go Xjj 





It can then be seen that this set of equations 2-11 may be vrritten in 
matrix form as! 





’ yo' 




go 


0 


0 0 


0 ® • 


® 0~ 




Xo 




7i 




gn 


go 


0 0 


0 


0 




Xj 




7z 




Sz 


Si 


go 0 


0 


0 




Xa 




ya 

o 

o 




S3 

o 

o 


Sz 

o 

c 


Si So 

O 0 

0 o 


0 ® • 

« 

o 


O o o 

o 




xa 








o 

Sn 


c 

Sn -1 


o o 

gn -2 gn-3 


o 

o o 


« 

• ®go 


L 


^n 


or 




y' 


1 = 


[“] 










where 


y; 


1 and 


i; 


1 are column matrices of 


order 


n 





and C® J is a square matrix of order n called the system trans- 

3 

fer matrix. 

For the case of a time-varying system one then obtains instead of equa- 
tion 2-7 and 2-8! 




g(nT,kT) x(kT) 



and the treuisfer matrix form for the system -G(s) would be! 

I 



(2-13) 



J 



9 



g(0,0) 0 



0 0 



0 




g(l,0) g(l,l) 0 

g(2,0) g(2,l) g(2,2) 

» o o 

» O O 



0 

0 



o 



0 



0 



( 2 - 14 ) 



g(n,0) g(n,l) g(n,2) g(n,3) » • » g(n,n) 

For example 1 in order to obtain the output response for a unit 
step input to a system vdiich is simply an integrator, we have to de- 
termine the X and G matrices. For a unit step input, the input at 
the sampling instants is always one or I 



X 




1 

1 

1 

» 

e 

1 



To determine the system matrix one must determine the weighting 
values at the sampling instants . If the period T is equal to one second 
then one simple method of determining the sequence is to find G(z) and 
to divide! 



G(z) = 



-1 -2 -3 

— 1 + z +z +z +°°® 



( 2 - 15 ) 



z - 1 

Another, more generally useful method is to determine the impulse 
response g(t) and find the values at the sampling instants. For an in= 
tegrator, G(s) = l/s or! 

g(t) = u(t) a unit step. 



Therefore, the system matrix is! 



[ 0 ] = 



10 0 0 
110 0 
1110 



0 0 
0 0 
0 0 



1 1 



( 2 - 16 ) 



10 



>3 




I 




10 0 1 

1 1 0 . . 1 



1 

2 



111 1 



3 

4 



(2-17) 



111 1 



2-4 The Evaluation of the ^stem Transfer ^fe.trix of a Transfer Function 
]ii order to determine the output response of a system, the system 
transfer matrix must be available. This is determined with the most 
facility and accuracy by determining the g(t) by taking the inverse trans- 
form of G(s), In order to evaluate the matrix values substitute t = nT, 
This method of evaluation of the matrix is illustrated in Appendix A 
and Table A-1 gives the values for some representive systems, 

Li many cases, a system has an undetermined G(s) or frequency res- 
ponse and the G(s) must be determined experimentally. In this case, it 
would be as convenient to determine the impulse response of the system 
directly and therefore the g^ values. For a large percentage of the 
sampled systems, a hold circuit filters the output of the sampler as 
shown in figure 2-3. In this case, the unit impulse is converted to a 
unit pulse of one sampling pe riod width. Hierefore, the z-transform 
equation for the output T (z ) is * 



Therefore, we are interested in determining the matrix 

experimentally the values of this matrix may be 
determined by exciting the system G(s) with a unit pulse of period T, 



Y (z) = Gj^G(z)X(z) 



where 




(2-18) 
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This method has been verified experimentally and yielded values within 
three percent of the expected elements of the natrix. This method of 
determination of the pulse response of the controlled system is easy to 
accomplish, and quite useful in investigating components with unknown 
transfer functions. 



;>-5 Response Between the Sampling Instants 

In any but the most well behaved system, the response between the 
sampling instants is of interest and must be determined. For this puT“ 
pose the modified z transform was introduced in the z domin, and an 
analogou? method must be determined for the time domain matrix, Ifriting 
equation (2-9) one has^ 

g„ = g(nT) = g(t)| (2-19) 

It = nT 

for the values at the sampling instants , Hierefore , if the values at 
half-v/ay between each sample are to be determined (m = 1/2 in the modi=> 
fied transform) one has” 



g (m = 1/2) - g(nT -T) = g(t)| 

It = nT iH -T 
2 



( 2 - 20 ) 



The symbol m was chosen to be cansistent with the modified z-transfonn 
and is defined as the percentage of the period from the sample point, 

as shovm in figure 2-4, The index m can assume a value 0 to 1 and is 

th. "th 

at the n sample when m = 0 and the n-1 sample when m = 1, 

For example, if the output between the sampling instants is required 

for the system shown in figure 2-5, then one writes the equation for the 

output as' 



Y(m)] = [G(m)] x] 



( 2 - 21 ) 



The impulse response of the system is ! 
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x(t) 


Hold 




G(s) 


y(t) 




<^(=) 







Figure 2-3 o Open Loop Sampled System V.fith Hold 




Figure 2-4, The Intersample Response 
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g(t) = e"^ 

Then, the values of the transfer matrix at the sampling instants 
(m = O) is° 



®n 

and the matrix is! 



-n 



[=] 



10 0 
-1 , 



-2 -1 , 
eel 



-n 



( 2 - 22 ) 



In order to determine the values at half-way between the sample 
points one must determine [o(m)J where m = 1/2 o Therefore 



=>n 



and 



[g(1/2)] 



-(n + 


1/2) 


“1/2 

e 


0 


“1.5 




e 


e 




e-1.5 



0 

0 



(2-23) 



Then, the output at the sampling Instants and at the mid-point of the 
sampling period ist 



C = 



10 0 


1 


1 


(.3679) 1 0 


1 


= 1.3679 


(.1353) (.3679) 1 

o 

• 


1 

o 

• 


1.5032 



(2-24) 



and 







".6065 


0 


0 


1 


0 

8 

1 


C(l/2) 




.2231 


.6065 


0 


1 


= .8296 






.0821 

• 


.2231 


.6065 

0 


1 


.9117 






o 


14 


0 




— 



(2-25) 
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The determination of these values is discussed further in Appendix A 
and values are given for the systems considered. 



2-6 The Formulation of the I-fetrix Equation for Closed Loop Systems 
Before considering the closed loop system, one must consider the 
two block open loop system as shown in figure 2-6 and examine the mat- 
rix algebra. 

Then, the necessary equations are' 

nj = [g,] x] and y]=[g 2 ]b] (2-26) 

Therefore, xj =[Ga] [g,] x] (2-27) 

and, in general, since matrix multiplication is not commutative it is 
incorrect to write the transfer matrices in the reverse order , that is ° 



y] = [g,][g,] x]/[g,] [g,] xJ 



( 2 - 28 ) 



Now, for a closed loop control system as shown in figure 2-7 one has! 
e(t) = r(t) - (c)t and (t) = r^ (t) - c* (t) 



In matrix form one obtains J 

e] - rJ- cJ and c] = [Gj] [g^J eJ 



(2-29) 



Therefore, one has! 



Therefore 



e] = r] - [G 2 ][g,] e] and e] = |l] ♦ [g,][g,]| ^ e] (2-30) 

<=] = [°<] {[ l] * [°a] [0<l] * "] (2-31) 



where [^l] = identity matrix = 



1 0 0 « • 0 

0 10 

0 0 1 •• 0 



« • 

a • 

0 . 



0 1 



and [aJ ^ = Inverse of matrix A. Thus, the solution for the output 

response involves matrix multiplication, addition, and inversion. 



Fortunately, one finds the inversion process is simplified by the fact 
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Figure 2-5. First Order Sampled System 



Y(s) 



:(t) 




y(t) 



Figure 2-6, Two Block Open Loop System 
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that the matrices are all lower trangular matrices; that is, all the 
elements above the diagonal are zero. Also, one finds equation 2-31 
may be written as * 

c]=[k]r] where [k] = [a]|i]+ [a (2-32) 

where [^] ” [^ 2 ] product of the transmission matrices 

in the forward path. 

Post multiplying equation 2-32 by ] ■•■ [A. ] one obtains ; 

[!']{[!] ^ [*]]= [a] (2-33) 

Adding “ L^] “ [*^] both sides, then; 

[K]{[l].[A]]=[Al*[l]- [l]=[[l]*[Al|-[l] 
Therefore, postmultiplying by the inverse of[l] + [A],one obtains ; 

[k]= [i]-{[i]*[a]|-1 i2-k) 

which substitutes subtraction for more difficult multiplication neces=> 
sary in equation (2-32). 

2-7 The Solution for the Response of a Simple System 

In order to evaluate the response of a closed-loop system, the in- 
verse of the matrix [l]+[A]must be determined. Since a physical system 
always is represented ly a lower triangular matrix, one powerful method 
for inversion is given by Frazer, Duncan, and Collar and presented in 
Appendix B.^ Consider a simple approximate first order transfer func- 
tion system where the sampling period is one second and a» 1, so that 
a neglible delay will be introduced, but there will be no output at the 
n = 0 sampling instant. The system with a step input is shown in 
figure 2-8. 

The system matrix is determined in Appendix A and given in Table 
A-1 and rewritten here; 
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[“] 



0 0 0 
10 0 



then 



[ll *[01 = 



1 

1 

1 



1 0 
1 : 
1 1 



0 0 
1 0 
1 1 



and 



{[II » ro}] 



-1 



1 

-1 

0 



0 

1 

-1 



Therefore , 



[I] = 



0 0 
+1 0 
0 1 
0 0 



0 0 

0 0 

1 ,, 0 

e 

0 ,. 1 

0 

0 

0 

1 . . 



so that 



o] = 



0 


0 


0 


1 


0 


1 


0 


0 


1 


1 


0 


1 


0 


1 


= 1 


0 


0 


1 . . 


1 


1 


[. 




• 

e 

o 

cJ 


e 

• 

1 


1 

0 

e 

1 ^ 



(2-35) 



( 2 - 36 ) 



(2-37) 



( 2 - 38 ) 



The values between the sampling instants after the first period 
would be the same as at the sampling instants, since 
[G(m = 0)] = [g(b)] for n > 1 . 
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II !; 





Figure 2-7, Closed Loop Sampled System 




Figure 2-8. Closed Loop System with A Step liiput 



19 



CHAPTER 3 



TflE MilTHEMATICAL THEORY OF INFINITE KATRICES 
3-1 The I'ktheraatical Theory 

It was shown in chapter two that the output time response could 
be obtained with the use of a matrix formulation. Furthermore, these 
matrices consist of elements whose values are those of a time response 
evaluated at discrete instants. For an open loop sampled-data system 
the matrix equation 2-12 will be rewritten here’ 

y] = [g]x] (3-1) 

The column matrices yJ and X | contain the values of the discrete 
response and input respectively. If the response is to be determined 
for all time, then the order of the column matrix, and the square sys- 
tem matrix, is infinite. Of course, even if one does not need to evalu- 
ate the infinite number of response values, the system matrix can be 
considered as an infinite matrix J that is, of large order. The content 
of this chapter is a discussion of infinite matrices and infinite se- 
quences, and the calculation of the time response of a ssimpled-data 
system utilizing a matrix formulation in the time domain. 

The system matrix G is an infinite matrix since 

[a] = A = (a^j.) (i,j - 1,2,3, 0 • . . . oo ) (3-2) 

that is, the matrix A is an array of elements of infinite order. Hence- 
forth, in this chapter, let [a] = A and the column matrices xj = x and 

"I 5 

Y j = y for convenience in notation, Cooke discusses, at length, the 

characteristics of infinite matrices. Only the characteristics of 

immediate importance shall be mentioned here. In general, the theory 

of infinite matrices' is connected with mathematical analysis and the 

theory of functions. 
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The operations of interest for the infinite matrices are° 
A + B = (a^ j , ) 



OQ 



AB 



= 2 



k=l 



^ik^kj 



(3-3) 



M = . 



where X is a scalar. 

The input and output time responses are expressed as an infinite 
sequence of discrete values. These discrete values constitute the ele- 
ments of the column vectors x and y, A vector x may be considered a 

5 

sequence space. A definition of sequence space is, 

A set S of sequences is called a sequence space when it contains 
the origin, and is such that, for every x and y in S and for every 
(complex) scalar c, x+y and cx are in S. 

The sequence space of interest is called a» the space of all se- 
quences. Then, the matrix equation 
y = Ax 
^n 

or y^ 



(3-4) 



a 



a , X, 
nk k 



k=0 



is a linear transformation of a on itself, that is, the system with 
matrix A transforms the input sequence into another sequence, the out- 
put sequence . When 

^nk ~ ^n-k (3-5) 

the system is time invariant as previously discussed in section 2-3, 

The infinite, positive time, ma. trices for physical systems are al- 
ways lower triangular mtrices (L.T.M. ), that is ; 
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a 



= 0 when j > i 



•\lsoj, for time invariant systems the elements are equal along the upper 
left to lower right diagonal „ This matrix is called a diagonally in- 
variant matrix (DoI.M,) and exists when" 



It is important to note, that a matrix multiplication is non-coramu- 
tative for the class of systems which are time- varying (T„V.)o These 
time-varying systems, for the operations defined in equation 3“3s ^ire 
characterized by a non-Abelian (non-commutative) algebra. The time in- 
variant systems may be called a sub-class (T.I.) of the large class T.V. 
and are characterized by a Abelian algebra. That is, for time invariant 
systems, matrix multiplication is commutative. Thus for linear constant 
coefficient systems the order of the natrix multiplication may be re- 
versed, while for nonlinear or time varying systems they may not. 

Inversion of the infinite matrix A is often necessary and it is 
shown in reference 3» that if there is no i where a^^ = 0, and A is non- 
singular, then a unique inverse of A exists. The evaluation of the 
inverse of the L.T.M. is discussed in Appendix B. 

3-2 Convergence and Stability in the Sequence Space. 

In an automatic control system it is of great importance to deter- 
mine if the overall closed-loop system is stable and therefore the out- 
put time response is bounded. The output response is a sequence in the 
sequence space and it must be determined if the sequence is bounded . 
Therefore, if a step function test signal is applied to the system it 
is necessary to determine if the output sequence converges to a final 
value. If the system response diverges, the system is considered unstable. 




for every i,j 

J 



(3-6) 
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For a closed loop control system, the output response may be written 
y = Ax (3“7) 

where A is the closed loop system natrix and 

A = G(I + G)“^ (3“8) 

The G is an open loop system matrix as shown in figure 3-1 

T 



X(s) 




System 

G(s) 



I(s) 



Figure 3-1 « Closed Loop Sampled Data System 
Then, given an input signal of a convergent nature, that is, of a 
bounded nature, the problem is to determine if the output sequence y is 
of a bounded nature, thematically, the problem is, given a convergent 
sequence x, under the transfonnation A, does a convergent sequence y re« 
suit ? A theorem of fundamental importance concerning this problem is 
that of Ko j ima-Schur . ^ 

Kojiraa-Schur theorem i the necessary and sufficient condi- 



tions that 



Ax or y = 
•'n 



CO 



k=l 






(3=9) 



should tend to a finite limit as n-»oo whenever Xj^ is convergent are that 



CO 



(a) 




nk 



^ M for every n- 



k=l 



(b) lim a^j^ - for every fixed k, 



n-^co 



( 3 = 10 ) 



(3=11) 
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(c) 




n->eo k=l 



(3-12) 



Moreover, if lim Xj^ = x then the final value is 

k-^oo 



CO 



(d) y = lim 



n 



X 




(x, 
k k 



X) 



(3-13) 



n -Hs* eo 



k=l 



A matrix satisfying a,b,c is called a k-matrix with and « for its 
characteristic numbers. For, the special condition of = 0 in condi= 
tion (b), one defines the matrix A as a matrix. These are the condi- 
tions for convergence and clearly define the stability of a control sys“ 
tem in the time domain. By examination of the impulse response of stable 
systems a physical understanding of these conditions results. Consider, 
for an example , a time invariant open loop system with a transfer function 

G(s) = and T = 1 as discussed in section 2-3» 



Then a^j^ = ®'’n-k ^m ~ I*®® where m = n - k. It is obvious that con- 
dition (b) yields = 1,00. Furthermore, it can be seen that this 
system does not satisfy conditions (a) and (c). That is, examining condi= 



tion (a) one finds for n 

Os 



that 



ZL hnk| 

K=1 



OD 



This result is as expected since the response 



is diverging as was determined in section 2-3» equation 2-17, Therefor* 
condition (a) requires that the series j®^nk| order for 

the system to be stable. For a closed loop system of figure 3-19 it is 
important at this point to recall that A = G o} , Therefore, A 
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may satisfy the condition (a) while G may not. 

Theorem" If A is the matrix for a time-varying system (T,V,), the sys- 
tem is stable if and only if it is a matrix. 

Corollary" If A is the matrix for a time invariant system, then it is 
a D,I.M. and A = , and A is stable if and only if A is a 

T matrix. 

OC 

For the large class of time invariant systems, the corollary implies 

that for stability it is necessary for 
CO 

< CO where n = i - j (3-14^ 

n=0 

For stable T,I. systems, one finds that = 0, that is the matrix 
is a matrix. For condition (c) one obtains « = 1,00 for a type I 
servo system. Then, by the use of relation (d) one obtains, as expected, 
the final value as’ 

y = « X = X (3“15) 

The transformation on a, accomplished by A as a T^ matrix is called 
a regular transformation. The property of importance is that as in 
equation 3-15 » the T^ matrices have the property of consistency. That 
is, every convergent sequence is transformed by such a matrix into another 
convergent sequence with the same limit when n approaches infinity. In 
addition, T^ or k matrices will frequently transform divergent sequences 

5 

into convergent sequences , On this point it may be stated * 

Corresponding to each unbounded divergent sequence { | each 

bounded sequence { y^^ ^ , there is a general (square ) T^ matrix which 

carries \ S I into i Y )• 

In this chapter, the mathematical basis for the infinite matrices 
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1 



has been discussed, and the operations defined. Furthermore, the theory 
of stability and convergence in the sequence space has been discussed , 

It is not to detract from this, that usually the conditions a), b), c), 
d) are investigated simultaneously with the evaluation of the actual 
output sequence for a given system. 
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CHAPTER 4 



ANALYSIS OF CLOSED LOOP CONTROL SYSTEMS 
4-1 A Ifethod of Evaluating the Response of a Closed Loop System With- 
out Inversion of Ife. trices 

Consider the simple, error sampled, closed loop system as discussed 
in chapter two and shown in figure 4-1. 




c(t) 



Figure 4-1, Error Sampled System 

It was sho\ifn that the output response time sequence may be written as° 

°]=M{[i] r] = +[g]|"^ k] (4-1) 

'Rierefore, the evaluation of the output response involves the inversion 
oftl] ♦ [O] as discussed in Appendix B, 

However, there is a simple method of evaluating the output response 
which avoids the inversion of a matrix. The output response column ma- 
trix may be written as " 

C]=[g]e] (4-2) 

where eJ = the error sequence in time . 

The error sequence may be written as 

e] = R] - C] (4-3) 

It can be seen that the error matrix can be easily evaluated by subtrac- 
tion. Then the error matrix multiplied by the system matrix will yield 

27 




the output response. Of course, the error mtrix cannot be evaluated 
until the output sample values are available. However, it nay be seen 
that the output sequence may be evaluated in a step by step procedure 
by examining the expanded form of equations 4-2 and 4-3 for the first 
three samplixig instants. For a linear, time invariant system, one has; 



(4-4) 



Co 




go 


0 


0 * 


Co 


Cl 




Si 


So 


0 








Sz 


Si 


go 


®2 

- 



and 


eo 


To 


Co 


O 

o 

o 




e, 

, 


= ^1 

. 


- c^ 
C2 


= (r, - c, ) 

(rj - Cj) 



(4-5) 



Expanding equation 4-4 with row by row matrix multiplication one obtains 
the following results. 

For the first sampling instant, the time origin, one has° 



and 



So ®o 



eo = 1*0 - Co 



i4-6) 



Therefore, solving for Bq one obtains! 



6o 



l~^go 



(4-7) 



However, for all physical systems, there can not be an output immediately J 
that is, there is always a time delay of small, but real magnitude. There* 



fore, go = 0 for systems of interest. Then one has 

and (4-8 ) 

Co = 0 

For the second sampling instant one has! 

c^ - gi6o + go6< - gi Bo = g,ro (4-9) 
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and 



e, = r, - Cl 

Since and Tq are known, then c, may be evaluated and then e^ nay be 
evaluated. Now, for the third sampling instant one obtains" 



S 2 ®o Si ®i So ®2 



S2®0 + Sl®1 



and 



( 4 » 10 ) 



®2 = ^2 



Therefore, it can be seen that one may proceed to evaluate the 
time response sample by sample in time , The equations such as 4-10 be™ 
come lengthy for larger n and it is simpler to use the matrix form. The 
matrix form may be written, with go = 0, as" 

®o 

c, = g^ u u Ir^ - c^ ) (4=11) 

®2 

o 
o 
o 

c 

n 

For a linear time invariant system, the system matrix is diagonally in™ 

variant. Therefore, multiplication of [ G] e] as a row-column multipli™ 

cation may be replaced by a column-column multiplication. This is possi- 

ble since, for example, the third row is identical to the first column 

rd> 

from the third element upwards. That is, one may write the 3 mul tipi id- 
eation as" 





~0 0 0 0 


To 




gi 0 0 


(r< - c^ ) 




g2 gl 0 

9 0 o 

o 0 c 


(^2 " ^ 2 ) 




o o 0 

Sn e^-1 Sn-2 


(r - c ) 
n n 





*®o ■ 




’®0 


®2 ~ 1^82 » Si 9 gflij 


®1 


= |g2» gl » 0 ] 


®1 




®2 

L ^ - 


- 


®2 



(4=12) 



or alternately as : 





”0 ' 




®0 


®2 = 


gl 




®1 








®2 



g2®0 + gl®1 



(4=13) 
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I-lultiplication is accomplished by multiplying the first element of 
E] by the last element of[G] of interest. 

Then, in general it is possible to write I 



[c] = 



*0 * 


eo 


0 


g1 


e, 


C1 


Sz 


62 


= °2 


63 




C3 


84 

0 


®4 

• 


°4 

0 


e 


0 


0 



(4-14) 



If a table of the form of equation 4-14 is established , a step by step 
evaluation procedure may be used as follows J 

1) Evaluate &q = Tq , and Cq = 0 

2) Then for c^ , multiply the G and E matrices by starting at g, 

in the G natrix. Then one obtains c, = gn e© , Wow evaluate 
e, - * c^i and fill in c, and 62 in the computation table „ 

3 ) Now evaluate C 2 by multiplying the G and E matrices by starting 

at gj. Then one obtains C 2 = Sz^o gi®i • Then evaluate 

®2 ~ ^2 ” °2 place the values in the table. 

4 ) Continue this procedure for each sample point of interest, 
typically until the system settles to a final value. 

An example will best illustrate this procedure. Let the system 
matrix of figure 4-1 be* 



0 ] = 



0 

.4 

.8 

1.00 

1.0 



0 

0 

.4 

.8 

1.0 



0 

0 

0 

.4 

.8 



0 

0 

0 

0 

.4 



0 

0 

0 

0 

0 . . 



( 4 - 15 ) 



Assume a step input signal =[l,l,l, . . . . ij , Then, es 
tablishing a table as in equat Ion 14 one calculates each value of c 



n 
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and step ly step as outlined. This method is illustrated for the 
first two calculations by use of dashed lines to indicate the flow of 
the calculations as follows ^ 

G E C 



c] = 



0 


A' ' “ 1»*^ 




Co = 0 


.4^' 


5, > ^ » 

e, = r, - c, = 1 - .40 = .60 


- 


c, = .40 




V-_ ^ 




✓ 


.8 






C 2 


1.0 


63 




C 3 


fcs C-3 









(4~16) 



Ttien, the next step is illustrated schematically ast 
G E C 



c] 



e, = 1.0 







Cq = 0 



o8* 

loO 



_ _ _ _ .i« 



1 - c, = - .040 1 



o40 

1.040 



- w: 



(4-17) 



e3 J L 

If this procedure is carried on step by step, one obtains for the 

response at the first six sample points *, 

G E C 



c] 



0 


1 


0 


.4 


.6 


,4 


.8 


-.040 


= 1.040 


1.0 


-.464 


1.464 


1.0 


-.3824 


1.3824 


1.0 


-.0358 , 


1.0358 



(4-18) 



To evaluate the Intersample response one proceeds as -discussed in 
section 2-5. Rewriting equation 2-21 one has* 

C(m)] =[G(m)j e] (4-19) 

If the G(m) matrix for the midpoint of the sampling period is founc 
to be! 
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1 



(4-20) 




,20 


0 


0 


0 


0 


,60 


.20 


0 


0 


0 


.90 


.60 


.20 


0 


0 


1,00 


.90 


.60 


.20 


0 



"nien one may calculate the intersampling response as° 



0 ( 1 / 2 )] 



G(m) 


s 


C 


p=> ■- 1 

.20 


• a: 

1 


,200 


,60 


.600 


.720 


.90 


-.040 


= 1.252 


1,00 


-.464 


1,423 




-.3824 


1.209 



(4-21) 



It is valuable to note at this point, that the calculations in the 



case considered in this section were simplified ty the unity feedback 
condition which yields e] ^ rJ ~ Sirigle loop systems with other 
than unity feedback shall be considered in the next section. 



4-2 Evaluation of the Response of A Closed Loop System with Other than 
Ihity Feedback, 

In this section the closed-loop system shall be considered with 
other than unity feedback as shown for one case in figure 4-2, The 
samplers are synchronized and with the same sampling period. Then the 
matrix equations my be written as " 

c] = [g] eJ ( 4-22) 

e] - r]» b] where b| ^[h]cJ 
Therefore eJ = rJ - [h] c] (4-23) 

Now, to evaluate the output response one uses equations 4-22 and 
4-23 and evaluates the and e^ step by step as outlined in the previous 

section. The only difference is in this case, in order to evaluate E^ , 
one must evaluate and then substract b^, the sample value, from 

r^ to obtain the error at the particular saTfiple, It must be pointed 
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Figure 4-2 . Two Sampler Feedback Ssrsten 




Figure 4-3 » Error Saiopled System 
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out that the solution in this case depended on the fact that the trans- 
fer blocks G(s) and H(s) were separated by samplers. 

Consider the system shown in figure 4“3 which does not have samplers 
separating both transfer blocks. Then on* may write" 

c| =[gJe| 

and e \ - R j « B] 

But, noting that B(s) - H ( 3 ) 0 ( 3 ) = H(s)G(s) E (s), then, 

k1 = r] - Hc] 

= h] - [ho] e] 

Therefore, it can be seen that e 1 is not readily available in equation 
4-26 since Hc] is not usually available. 

Then, rearranging equation 4-27, one obtains^ 






(4-24) 

(4-25) 

(4-26) 

(4-27) 



f [ll - H } e] = h] 



cr e] =([i] * [l 



HG 



R 



( 4 - 28 ) 

(4-29) 



Therefore, in order to evaluate the output response for this system, 
inversion of a mtrix cannot be avoided. 

Finally, if one was analyzing a system such as shown in figure 4-4 9 
the equations of interest are" 



E 



R 



[h] c] 



ana c] = [g,1 [g,| e] 



(4-30) 

(4-31) 

If the response at an inteneediate point in tiie system is desired , such 
as the output of G, (s), it is readily available by writing the following 
equation" 

Ea] - [g/] e] (4-32) 

The response at intermediate points In a system is of great importance 
in nonlinear systems which ifill tie discussed in section 4-5, 
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Section 4-3 Analysis of Multiloop Control Systems „ 

The Introduction of more than one feedback loop is necessary or 
inherently present in many control systems. The investigation of a 
multiloop sampled-data control system is complicated by the presence 
of samplers in some loops, and the absence of samplers in others. There- 
fore, all the feedback signals are not of the same form throughout the 
system. Consider at first a two loop system as shown in figure 4-59 
which has samplers separating all the transfer blocks. 

One my write a set of simultaneous s;impled-data equations, where 
the starred notation indicates a sampled signal or transform, as follows t 
(s) - n (s) - h"' (s)c'^(s) - n (s) 
m'^Cs) = G^(s)e"' ( s) (4-33) 

C*(s) = G^(s)M*(s) = G*(s)g;(s)e"(s) 

Solving these equations simultaneously, one obtains for the sampled- 
output ° 



C (s) 



G^(s)G^(s)R*(s) 

1 + G*(s) + g;(s)G^(s)H*(s) 



(4-34) 



The 2 -trans formed equation then follows as° 

C(z) ■■ 



G^(z)G^(z )R(z) 



1 G^(z) + G, (z)G (z)H(z) 

a Da 



(4-35) 



One nay then write the matrix equations directly from equations 4-34 
and 4-359 or alternatively one may derive the matrix equations directly 
from the system signals. In either case, one obtains J 

o] = |[ ll * [ * [Obl [\1 M [oJ R] (4-37) 

«d 0(m)] = ([l] > [gJ * [oj [ qJ[h]I ^ [Gj,(d)] [g^(d)] r] (4-38) 
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Figure 4-5. Multiloop Sampled System 



36 



If it was desired to avoid the inversion necessary in equations 4-37 and 
4-38, one may write a set of simultaneous equations which provide a means 
of step by step evaluation of intermediate signals. In this case one 
would writer 



e1 = r 1 - Ml - [hJ c] 


(4-39) 


»]= KM 


(4-40) 


=] = K1 «] 


(4-41) 



In order to evaluate sj , m], and a step by step procedure is 
followed similar to that of the single loop method. First, assuming 
no immediate o'Xput for either or , one may i^rite^ 

Bo " 0 

c© - 0 (4-42) 

and 6o = ro 

Then, for the next sample one obtains* 
a *0 “ a ^0 

vdiere g^g. “ first impulse response value of the G block and 

8i 

= the n^^ impulse response value of the G^ block. 

Then, continuing, one obtains* 

= g^b ^ gob = gib ®o 0 (4-43) 

since ggb ~ 0 and nio = 0. 

Therefore, one uses equations 4-39, 4-40, 4-41 in that order, step by 
step. For the next sample one obtains* 

e, = r, - m, « (hgC^ h, Cq ) = r, - g, g Tq - 0 

(4—44) 

since Cg ~ c^ - 0 

and ffl 2 = gaa eg + g,g e, (4-45) 

and finally 
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= ggb -^gib^ ”gnb“^ since mo - 0 

As is the usual case, it is actually easier and more methodical to use 
the matrix form throughout the calculations , Hieref ore , one vrill esta- 
blish a table similar to that of section 4”1 using equations 4-39 » 4-40, 
4-41. The flow of calculation is from e1 to to c] and then back 
through the feedback loops again to eJ , It is obvious, that actually 
one is simply following the signal around the closed loops. 



Finally, consider a multiloop system which possesses only one samp- 
ler in the control loops. The simple sampler typically is placed in the 
error channel as is shown in figure 4-6, 




Flgme 4“6. A Multiloop Sarapled Sj'^stem 

One may write the equations for the signals as follows ° 

E = R _ M - B 
M = G E^ 

^ ^ (4-46) 

C - G, M = G, G E 
D Da 

B = HC 

Therefore, one obtains for the closed loop sampled error and output' 

1 + g'^ 

a b a 






■s? 

E 



(4-47) 



and 



X 

C s 



(4-4B) 






1 + G 



R 

HG, G 

b a 






One may then vrrite the matrix equations directly as* 

=[[l]"[°al »] 

= [°a“bl[M *[°al ^ 

The equation for the intersample response is then* 

C(m)] = [0^Qb(n.)l[[l] *[gJ *[hG^0J|-^ r] 



(4-49) 

(4-50) 



(4-51) 



Thei^fore, whenever one desires to calculate the intersample response 

and the inversion of the mtrix I + G + HG, G has been already accom- 

a D a 

plished, one simply evaluates the G G, (m) matrix and carries out the 

a u 

multiplication. Therefore, multiloop systems with one sampler or samp- 
lers separating all transfer blocks may be equally treated by the use 
of time domain infinite matrices. The response at any sanqjler location 
is readily available and is usually of interest in the investigation of 
multiloop systems, particularly with nonlinearities present. The intro- 
duction of nonlinearities into a control loop is treated in section 4-5 o 



Section 4-4. Analysis of Time-varying Control Systems 

The analysis of time-varying sampled-data control systems may be 
accomplished using the tine domain matrix method. A control system 
may have as one of the transfer function blocks, a component whose 
parameters are changing with time , This effect is present with high 
altitude jet aircraft, where the dynamic characteristics of the air- 
craft change with altitude , and therefore time , 

Consider the simple open loop system shown in figure 4-7, The 
impulse response of the time varying component is changing with time . 
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Tine-Varyixig 
Component G 



c(t) 



Figure 4-7 , Sampled Open Loop System 



Therefore the matrix equation for the output may be written as in section 
2-3, equation 2-13 and rewritten here" 

c] = [G(n,k)] r] (4-52) 



where 



[G(n,k)] 



g(0,0) 

g(l,0) 

g(2,0) 



0 

g(l,l) 

g(2,l) 



0 

0 

g(2,2) 



e e o 



o o o 



In order to substitute numbers in the G(n,k) matrix, the time variation 
of the system G must be known or determined. As an illustration, consi- 






der an abrupt change of 

G(s) “ from a = 1,0 to a = 2,0 at the third 

S 'r S. ) 

sampling instant (k ~ 2) where T = 1 second. Then, using the values of 



from table A-1 one has for the system matrix^ 



[o(n,k)] 



0 


0 


0 


0 


0 


Q e e • 


3679 


0 


0 


0 


0 


0 


7675 


.3679 


0 


0 


0 


0 


9145 


.7675 


.2838 


0 


0 


0 


9685 


.9145 


.4707 


.2838 


0 


Q o o « 


9884 


.9685 


.4960 


.4707 


,2838 


0 
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Then, of course if this component had feedback introduced as shown 
in figure one would have the following closed loop equation" 



c] = [o(n,k)l |fl] . r] 



(4-53) 




c(t) 



Figure 4-F, Time Varying Closed I^oop ^/stem 
In order to obtain the intersainple response, the equation for the output 
response may be written, 

= ("’] ♦[0(n,k)]]“^ R] (4-54) 

Analysis of a feedback system with more than one tixte varying element 
follows the same approach as for the time invariant systems, Ite matrix 
equations are found to be the same as for time invariant systems with 
the time variation of a transfer function only affecting the system ma= 
trix itself. Therefore, for the system shown in figiare 4=-9, the follow- 
ing equation is obtained^ 

c] = [G(n,k)] |[l] + [G(n,k)l [H(n,k)]|“^ r] (4-55) 

It can be seen that it is not necessary for samplers to separate the 
time varying component from all other transfer blocks. Therefore, if 
H(s) was time invariant, it would not be necessary to have a sampler 
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Time Varying 






H(s) 





Figure 4“9o System with Two Time Varying Components 
between G(s) and H(s) for the use of the time domain matrices. For 
the system shown in figure 4~10, one obtains the equation, 

C] s [G(n,k)]^[l] > [GH(n,k)]| r] (4-56) 




Figure 4-10, Time Vaiying System with A Feedback Component 

A control system, where one of the components is a time=varying 
amplifier is worthy of consideration. Consider the system as shown in 
figure 4-11, 

The gain a(t; is changing with time and therefore has a different 
magnitude at each sampling Instant, The equation for the output may 
be written^ 
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Figure 4-11 . Open Loop Time Varying System 
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Then, the output at the third sampling instant is as expected; 
Cg = aogaro + a^g^r, 

It can be seen that equation 4-58 may be written as follows" 



V 


”0000^ 




0 

0 

0 

0 


Oj =[g]Wr] = 


• 

• 

0 

0 

0 




0 a, 0 0 




ga gi 0 0 




0 0 aj 0 




g 3 ga gi 0 




0 0 0 ajj 






(4-57) 



(4-58) 



(4-59) 



(4-60) 



The operation of multiplication of the A and G matrices is not commu- 
tative as was previously discussed „ That this is so is obvious from 
the form of the matrices. Now, consider an open system identical to 
that of figure 4-11 except that a hold circuit immediately follows the 
sampler, Tben, if the amplifier changes gain during the sample period, 
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this must t)« accounted for in the system matrix G, which is a function 
of m. In most cases it is reasonable to make the simplifying assump= 
tion that the gain changes only at the sampling instants. In other words 
the assumption is that if the gain is changing continuously with time, 
the time constant of this change is much greater than the sampling period 
Therefore, one is approximating the gain by a staircase function as 
shown in figure 4-12, 




t^ 



Figure 4-12, A Time Varying Gain 

Then, on this basis, one may analyze closed loop systems with a time 
varying gain. For the system shown in figure 4-13, one may write the 
equations for the system signals as’ 

e] = r] - [h] c] (4-61) 

and C] s [g] [a] e] (4-62) 

Equations 4-61 and 4-62 may be solved by the step by step procedure 
previously outlined for the time-invariant systems. Using this step by 
step method, the gain change at each interval is clearly displayed to 
the investigator. Alternately, one may solve for the output response 
and error response by means of inversion of the matrix in the following 
equations ’ 
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Figure 4-13. System with A Time-Varying Gain 

e] = {[i]*[»][<=] W] RJ (4-63) 

cj = [g][a]{[i>[h][g1[aJ|-1 R] (4^) 

If there was no sampler between the output anc3 the fee(3back block 
H(s), then one obtains the equation for the error sequence in the same 
manner as carried out in section 4-2. Therefore, one obtains for the 
error sequence. 

e] =[[i] * C«o] R] (4-65) 

As an example, consider a system where H(s) =1 as shown in figure 

-it 

4-13. Let R(s) = l/s, and the gain function be a(t) = e ^ , with a 
sampling period of one second. In this the time constant of the gain 
change is only double the sampling period, but the results are instruc- 
tive with this marked gain change as shown in figure 4-14. 

Use equations 4-61 and 4-62 where f hJ = [l] and one may write* 

e] = r] - C] (4-66) 

c] =[a][A>]= We] («’) 

where p] = j^Aj e] the amplified magnitude of the error pulse 

and 
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Figur« 4-14. Time Varying Gain Example 
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Then, one obtains^ 
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• 



(4-68) 



The availability of the amplified magnitude of the error pulse as p] 
is often useful in the investigation of nonlinear systems which shall be 
discussed in the next section. 
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4-5. Analysis of Control Systems with Nonlinear Components. 

The analysis of sampled-data control systems with nonlinear com- 
ponents may be accomplished by the use of time-domain infinite matrices . 
The methods most commonly used for continuous nonlinear systems are the 
describing function suid the phase-plane. These methods are somewhat 
limited in their application to sampled-data control systems. Consider 
a single loop nonlinear system as shoim in figure 4-15. 




Figure 4-15. Single Loop Nonlinear System 
Then, it can be seen that the nonlinear component has replaced the varia- 
ble gain amplifier of figure 4-13. The nonlinearities considered are 
those sensitive to the magnitude of the signal input J that is, they 
have a nonlinear amplitude response such as shown in figure 4-16 for a 
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Therefore, as the input voltage varies in magnitude , the gain of the 
amplifier changes. Since the input to the nonlinear element is the 
held eiTor magnitude, this magnitude determines the gain. However, the 
error is evaluated directly in the time domain and therefore a gain cal- 
culated for a calculated error input is also known as the gain at the 
specific sample Instant, Therefore, the nonlinear element may be ^ 
treated as a time-varying gain. This statement is true for hysteresis, 
a relay servo, and other single nonlinearities following a sampler. 

Then the output of the nonlinearity is the sequence , One may write 
for the output response and error response^ 

(4-69) 
(4-70) 

0 • • • 

0 

0^*9 

• o • 
e 

The values of u^ depend upon the error magnitude and are determined in 
the step by step solution. As an example, consider a saturating ampli- 
fier with a characteristic curve as shown in figure 4-17, T«<hich is the 
nonlinear component of figure 4-15, If the input signal is a unit step, 
then the amplifier is expected to saturate. This amplifier has a maxi- 
mum gain of two throughout the linear region. Then, using equation 
4-69, at the first ssuapling instant (n = 0)eo“l,0 and therefore 
u© ~ 1«0» and p© = 1,0, a gain of one. The G matrix for this example is" 



e] = r] - c] 
0] = [g] [0] E] 

where 



u© 0 0 

0 U) 0 
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Figure 4“17, Saturating Amplifier Characteristic Example 
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Therefore, c, = g^eQ = ,40 and e, = .60, This magnitude continues 
to saturate the amplifier and p, =1.0, and u, = 1,667. Continuing in 
this manner, one obtains the following matrix solution t 

c] = [o] [U] e] = [g] pJ (4-71) 
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The analysis of systems with a nonlinear component which is also 
a storage device follows along in the same manner. For example, if the 
nonlinearity was hysteresis in a magnetic amplifier, then the gain is 
dependent upon the past value of input voltage to determine which side 
of the hysteresis loop is applicable to the present signal. This can 
be seen clearly from figure 4-18 which shows a simple hysteresis loop. 

If the input was , there are tvro possible output magnitudes depending 
upon the magnitude of the input of the previous sample. If the previous 
sample had a magnitude of Eq, then the output magnitude will be E^, 




Furthermore, any nonlinearity may be treated in this manner since the 
relation 

P] = [bJe] (4-72) 

applies to ar^ nonlinearities whose characteristics are knovm or can be 
approximated , 

The analysis of systems with a nonlinear component which does not 
have a sampler immediately preceding it cannot be treated by this method. 
However, as an approximation a fictitious sampler may be placed before 
the nonlinearity, Ihe accuracy of this approximation depends upon the 
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sampling rate and the time constants of the system. This approximation 
shall be discussed further in^ the next section , 

If there is more than one nonlinearity present in the system and 
samplers appear before every nonlinearity, then the investigation of the 
system by time domain matrices is entirely possible. In fact, for a 
system with many nonlinearities, analysis is possible with the introduc- 
tion of san^jlers with a high sampling rate, wherever a nonlinearity 
exists. The necessary sampling rate shall be discussed in the next sec- 
tion. 

Furthermore, the signal sequence magnitudes are available at inter- 
mediate points throughout the system and can aid an investigator in the 
analysis of the component requirements, Consider for example the system 
of figure 4”19 which has two nonlinearities , Then the equations for the 




Figure 4-19. System with Two Nonlinear Components 
error, the output sequence of G, , and the output sequence of G 2 is! 



e] = r] - c] 


(4-73) 


H] = [0,1 [0,] E] 


(4-74) 


<=] =[0j1Wm] 


(4-75) 



A phase plane portrait of the error and the derivative of the error 
often aids an investigator in the analysis of a nonlinear system. By 
the use of the backward difference formulas of numerical methods, the 
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derivative of the error may be obtained using the past values of the 
discrete error. The derivative of the error nay be calculated ly use 
of the present and one past value asl^ 



e» = 
n 



(«n - 



(4-76) 



where the inaccuracy is of the order of magnitude of the sampling period 
T 

one has 



If the present value and two past values of the error are to be used, 
6 



'A = ^ * 0(tM (4-77) 

Use of three values of error has reduced the inaccuracy of the approxi- 
mation to within the order of the sampling period squared , order for ' 
this calculation to be accurate, the sampling period must be short with 
respect to the time constants of the system. This requirement is not 

I 

restrictii^ since this condition is necessary for stability in a sampled- 
data control system. 

If it was useful to postulate the phase space and determine the 
derivatives of higher orders it is possible to use the following formu- 
las, the choice of formula depending on the accuracy required,^ 



n 

®n = 


1 

T2 


('n - 


2*n-l ^ Va' * 




(4-78) 


/// 

*n = 


1 

T3 


(*n- 


^n-1 * 3»„-2 - 'n-3' 


+ 0(T) 


(4-79) 


/! 

*n 


1 

T2 




- 5'n-l * -‘Va - 'n-s' 


+ 0(T2) 


(4-80) 



Formulas for higher derivatives using backward differences are readily 

7 

available , 

Writing the above formulas in a matrix form results in the follow- 
ing matrices with an accuracy of the order of the sample period. 
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(4-81) 
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(4-82) 



The matrices vfith an accuracy of the order of the square of the sample 
period are* 



=] 
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-4 
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-4 

-4 




(4-83) 




4-6, Analysis of Continuous Control Systems 

The analysis of continuous control systems can be accomplished with 
the use of time-domain infinite matrices. This is possible through the 
introduction of a fictitious (mathematical) sampler or samplers in the 
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continuous closed loop. The approach is basically that of the use of 
numerical analysis methods in the solution of a differential equation. 
Therefore, the accuracy of this approximation depends upon the sampling 
rate of the fictitious sampler. However, the fact that it is not neces- 
sary to calculate the closed-loop roots or to invert the closed-loop 
response equation in order to obtain the time response, is of great im- 
portance o 

Consider a continuous system to be investigated, which is a single 
loop system as shown in figure 4-20. 




Figure 4-20, Continuous Feedback Control System 

2 

The approximation introduced by the fictitious sampler will depend upon* 

1) The location of the fictitious sampler 

2) The form of the fictitious hold circuit 

3) The frequency of the sampling 

Since it is important for the sampler frequency to be many times greater 
than the highest frequency of the input signal to the sampler, the loca- 
tion of the sampler is usually chosen to be in the feedback loop before 
the H (0 as shown in figure 4-21, This location takes advaintage of the 
filtering of the input signal R(s) by the plant G(s). 



54 




Figure 4-21. Location of the Fictitious Sampler 

Ihe fictitious hold network is present in order to reconstruct the 
continuous signal from the sampled signal. In actual sampled-data sys- 
tems a zero-order hold is usually used for its practical realizability 
and for stability considerations. However, for a fictitious hold the 
possibility of straight line and parabolic approximtions should be con- 
sidered, Consider a straight line approximation hold which can be 



achieved by a triangular hold circuit which has a transfer function ^ 



,2 
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The approximation of a time function is shovm in figure 4-22, This hold 
is not i^ysically realizable , but is very useful for mathematical appro- 




Figure 4-22. Hie Approximation of A Time Function 
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ITie sampling frequency determines the accuracy of the mathematical 
approximation and must be balanced by the ajnount of calculation that will 
be acceptable to the investigator. As the sampling period approaches 
zero 9 the inaccuracy is approaching zero, but the number of calculations 
rapidly approaches infinity. Fortunately, it has been determined experi- 
mentally that there exists a reasonable sampling rate for closed loop 
systems which yields less than 5% error in approximation « For a type I, 
second order system, this approximation holds when the sampling frequency 
is ten times the magnitude of the pole of the plant G(s), As an example, 
for a plant with a transfer function J 

the sampling frequency should be 10 radians per second or the sampling 
period is approximately one-half second. The calculations of previous 
sections were carried out with a period of one second and a type’ one sys- 
tem, Therefore, the number of calculations necessary to determine the 
approximate response to a continuous-data system are not impractical, 

Ihe simplest method of determining the necessary sampling rate is trial 
and error. One calculates a few points on the output response using a 
trial saunpling rate, then recalculates these points using twice the 
sampling rate. If there is a neglible change in the results, then the 
former rate vrais sufficient for the desired accuracy. 

This method of approximation may be used on non-linear or time- 
varying systems and therefore all the previously discussed methods of 
analysis will apply. 
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CHAPTER 5 



INTRODUCTION TO THE DESIGN OF CLOSED LOOP CONTROL SYSTE^6 
5-1 . Litroduction 

The method of time-domain infinite matrices may be used success- 
fully for design of closed loop control systems. The design or develop- 
ment of control systems for a specific application is a practical pro- 
blem of great importance. The designer may readily apply the analysis 
methods of the previous chapter to the design problem. They are used 
to evaluate the performance of the system under the specified operating 
conditions. Also, with the aid of experience in calculating the response 
directly in the time-domain, the designer may determine what system 
parameters must be adjusted sind in what manner. 

Given the basic specifications of the system and the performance 
requirements, the designer often must determine a compromise between 
conflicting requirements. In the design of closed loop continuous sys- 
tems, many designers rely on the open and closed loop frequency response 
curves of the system to indicate the performance of the system. These 
frequency response techniques are not very useful in sampled-data sys- 
tems. Correlation of the time response of the system with the frequency 
response is not readily achieved since a transformation of the z varia- 
ble, into a new complex frequency variable w, is necessary to map the 
unit circle of the z-plane into the entire left half of the w plane. 
Therefore, the use of frequency response characteristics, such as the 
height and frequency of the resonant peak, and the bandwidth, is very 
limited. Design by use of the root locus method in the z-plane is 
limited since l) the design only considers the response at the sampling 
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instants , and the inter sampling response may be wholly unacceptable , 

2 ) Correlation theorems between the time domain and the z-plane are 
accurate only under certain conditions . 

The time domain performance indices such as rise time , settling 
time, over-shoot, and number of oscillations are readily applied to 
sampled-data systems using the time domain matrix method. Furthermore, 
this method allov/s one to design directly in the time domain and deter- 
mine the compensator necessary for the specific application. This method 
of design of a digital compensator shall be discussed in chapter 7, 

For the design of continuous systems , the designer has all the 
standard techniques at his disposal for the selection of a system adjust- 
ment or compensator. The methods of Bode, Nichol, and the use of the 
root locus all may be applied in the design of the continuous system. 
Then, the sampled data approximation may be introduced and the design 
evaluated directly in the time domain by means of the time domain matrix 
method. Therefore, all the background of previous determined theory and 
laethods will apply profitably. The designer may use these techniques 
to determine the type and location of the compensator, and the parameter 
values . 

The designer may determine the necessary compensator directly in 
the tine domain by use of the digital compensator design technique dis- 
cussed in Chapter 7. If the system is a continuous system approximated 
by a sampled-data system, the digital compensator may be converted to a 
continuous compensator by time domain network synthesis techniques , 

In addition to these design problems, one may consider the design 
of adaptive systems. Due to the great interest of these systems, they 
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1 



are reserved for discussion in Chapter 6, This chapter shall be con- 
cerned with the use of standard design techniques and the use of time 
domain matrices to design various typos of control systems such as non- 
linear and time varying systems. 

5-2. Tlie Design of A System Utilizing Time Varying Gain Compensation. 

Consider the design of a closed-loop, unity feedback sampled -data 
system as shown in figure 5-1. The basic design steps are. 




Figure 5-1. Unity Feedback Sampled-Data Control Slystem 

(1) Determine the response for the given system G(s) and compare 
with the required perfonnance specifications. 

(2) If the performnee is not satisfactory, adjust the gain or 
choose another plant if possible , Otherwise , introduce a compensating 
element or block in the closed-loop. Choose the compensating component 
on the basis of factors such as design criteria and design experience 
and the time response obtained for the uncompensated system, 

(3) EvaliMite the compensated system response and readjust the sys- 
tem parameters if necessary. 

In this section it shall be assumed that 1) and 2) have been accom- 
plished and it has been decided to attempt to compensate using a time- 
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varying amplifier in the forward transmission path as shown in figure 
5-2, The variation of the gain must be chosen to provide the desired 
output response performance. Changing the gain of the an5)lifier does 
not alter the system dynamics and therefore the achievable results are 
limited. The designer may alter the response, but it is only possible 
to compromise between desired performance indexes. For example, it is 
possible to reduce the rise time, but only with a resulting increasing 
maximum overshoot. 
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Figure 5-2, Ihe Compensated Control System 



An example will best illustrate the possibilities. Consider a sys- 



tem with a transfer function of 



G(s) = 



1 

s (s + 1 ) 



(5-1) 



and a zero order hold and sampling period of one second, Ihe input sig- 
nal for this example will be a unit step. Then the imcompensated output 
at the sampling instants is found to be I 
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Therefore, the uncompensated system has a rise time of two seconds and 
an overshoot of greater than 40^ occurring between the third and fourth 
sample. If the specifications call for a maximum overshoot of 30 and a 
rise time less than four seconds, the required response can be achieved 
with a time-varying gain. The designer learns, from the experience of 
calculating the response directly in the time domain, that the first two 
error samples largely determine the magnitude of the maximum overshoot. 
Therefore, the designer choses an amplifier with a gain of one-half at 
the first two sauries, and a gain of one thereafter. This amplifier can 
be practically realized 1:^ constructing an amplifier which switches to a 
gain of .5 when a step input is applied, and switches to a gain of 1,0 
after two seconds ( n = 2), This system will give a more desirable res- 
ponse than simply lowering the gain to one-half for all time. This fact 
shall be verified later in this section. Calculating the compensated 
response one obtains' 

0] = [0] [Aj E] = [o] p] (5-3) 

and the refold 
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The intersawplp response is found as follows* 



C(ni)] = [G(ro)] P] = 



.10653 


.5 


.0533 


.6166 


.4080 


.3517 


.8590 


.4662 


= .7307 


.9481 


.0581 


1.1817 


.9809 


-.2365 


1.2884 


.9930 


-.2733 


1.2137 



(5“6) 



It can be seen that the maxiraun overshoot is about 



29% 



for the compensated 



system and the rise time is 3.15 seconds. Ibe response of the uncompen- 
sated, compensated, and the system with a gain of one-half for all time 
are shown on figure 5-3. A comparison of the rise time and maximum over- 
shoot is given in table 5-1. 



System 


Rise Time 


I*bximum Overshoot 


(seconds) 


(Percent) 


Uncompensated (Gain =1) 


2.0 


45 


Time -varying Gain 


3.15 


29 


Gain = ,5 for all time 


3.70 


15 



Further discussion of the application of the time domain matrix method 
to the design of control systems is presented in the next section, 

I 

5-3. Application of the Time Domain Jh.trlx ^fethod to the Design of 
Various Types of Control Systems, 

In this section it is intended to show the great scope of possibi- 
lities of design of various types of control systems using the time-domain 
matrix method, 

a) Single Loop Linear Sampled-Data Control System 

When the preliminary Einalysis of a linear sampled data system reveals 

•r 

that the overall transient performance is inadequate, compensation techni- 
ques must be employed in order to improve the system performance. For 
a single loop system, the simplest and most direct step is to change the 
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system gain. Usually, however, this adjustment alone is not sufficient 
to satisfy the design requirements. Therefore, it is necessary to in- 
sert a compensating network in the system in order to achieve the desired 
response , This compensating network may be cascaded in the forward or 
feedback channel or inserted as a minor feedback or feedforward loop. 

The sampling process in the control loop complicates the choice of the 
location for the compensating netv/ork. For error sampled systems, it 
is found advantageous to operate on the sampled error by a cascade com- 
pensator as shown in figure 5-4. 




Figure 5-4. Continuous-Data Compensation 

This compensator operates on the sampled and held error and yields con- 
tinuous data information to the controlled system. Furthermore, for the 
output response one may write" 

c] = [Gh^c^l e] (5-7) 

where is the system matrix for G^ (s )G^ (s )G (s ) . Investigation 

of this form of compensation reveals that it is difficult to stabilize 
a sampled-data control system containing higher-order integration with 
the use of linear continuous-data networks. The stabilization aiid com- 
pensation of sampled-data systems by means of continuous cascade 
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compensation is further complicated by the calculations necessary in 
order to find the natrix for the overall system for every trial compen= 
sator. This is not a disadvantage peculiar to this method, but also 
results when using the z-transform. The use of a continuous-data net- 
work for compensation is not to be excluded , These networks are simple 
R-C networks coupled with amplifiers and are easy to realize „ Further- 
more, for simple systems, they are sufficient and therefore probably 
desirable. However, for more complex systems it is often necessary to 
evaluate many trials in order to arrive at a reasonable compensation 
design. 



Therefore, it often becomes desirable to use a sampled-data network 
as a cascade compensator as shown in figure 5 - 5 „ Then, for the output 
response one obtains ' 




Figure 5-5. The Digital Compensator in A Sampled-Data System 

'=] = [°h“] [®] f] (5-8) 

where = Qd] e] (5-9) 

Therefore, the sampled input to the controlled system is the output of 
the compensating device which is a result of a transformation of the 
error sequence. As was stated in Chapter 3, there is always a [d] which 
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will transform a diverging sequence e] into a stable, converging se~ 

quence P] . Furthermore, the design of [dJ is relatively easier tiian 

the design of a continuous data compensator since the matrix need 

be evaluated only once . "nie samj^led-data compensating network can be 

realized as a program of a digital computer or more simply as a data- 

processing network. It has been shown, that a digital processing unit 

g 

may be realized ty operational amplifiers and electronic samplers „ 
Therefore, a digital compensator may be part of a large computer program 
for missile control, or simply a controller for a DC motor. The realiza- 
bility and the synthesis of digital compensators shall be discussed 
further in Chapter 7. 

B) Miltiloop Sampled-Data Control Systems. 

The design of multiloop sampled-data control systems is more com- 
plex and difficult than design of single loop control system. The use 
of linear continuous data compensation networks has the same limitations 
as discussed in- the previous paragraph. Therefore, the sampled-data 
compensating network is used more often. A fundamental problem in the 
design of multiloop systems is the selection of the location of the com- 
pensators. This problem is difficult to solve in continuous systems as 
it is in sampled-data systems . One advantage of the time domain matrix 
method is the availability, in the calculations, of the response at 
intermediate points throughout the multiloops, wherever a sampler exists. 
Kiowledge of each response allows the designer to use this informa.tion 
to adjust the parameters of the compensator. If there is only one 
sampler present in the multiloop system, then this advantage is not pre- 
sent. If a pulsed-data network is used as the compensator, then often 
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feedforward control wil] aid in the elimination of disturbance inputs „ 
Consider figure 5-6 which shows a feedback control system with a distuT” 
bance input and feedforward control . Then the output due to the inmt 
signal C, and that due to the disturbance signal can be considered 
separately as^ 




Figure 5-6, Control System with A Disturbance Input 

cj = [g] p] (5=10) 

Cj ={[I] 

Where P] = [Dj] R] + [d,] r] 

Then, the output due to the reference may be written as, 

c] =[o]|[D,] *[d,]|[[i] .[d,] R] (5-12) 

Therefore, the output due to the disturbance may be minimized by means 
of , and the digital network Dj used to design the system for the out- 
put performance with respect to the signal input. Multiloop systems will 
be considered further in the next chapter, especially the conditional 
feedback system, 

C) Nonlinear Control Systems. 

The techniques used in analyzing nonlinear systems as presented in 
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the previous chapter are very useful in the design of nonlinear systems o 
The nonlinear component is considered as a time-varying amplifier in the 
system where the gain is dependent upon the input magnitude , The design 
of linear as well as nonlinear compensators for the nonlinear systems 
is possible „ Also it is usually an advantage to have available the sig=> 
nal magnitude at various intermediate points in the system. 

The first method of compensation to be considered would be the in- 
sertion of a linear compensating network at a point in the system located 
before the nonlinear element. If this proved to be unsatisfactory, then 
a pulsed data network could be inserted before the nonlinearity in the 
loop. Also, since the nonlinearity is considered as a magnitude sensi- 
tive device and treated as a time varying amplifier, an interesting possi- 
bility for compensation would be the use of a time-varying amplifier with 
essentially complementary characteristics to that of the nonlinear device. 
If the amplifier and the nonlinear element were cascaded, then the magni- 
tude of the signal into the nonlinear element could be kept within the 
linear region of the component. In the case where the nonlinear element 
has storage of energy, this could not be achieved by a time-varying am- 
plifier. Therefore, in most cases a pulsed-data network is used as a 
compensator and designed on the basis of performance criteria and eval'ua- 
tion. 

The design of a nonlinear system usually involves some trial and 
error steps in order to evaluate effective design changes. The design 
of a relay servo, for example, may require a number of trials in order 
to arrive at the proper output voltage on the relay and acceptable values 
of dead-zone and hysteresis . The use of a nonlinear element as a com- 
pensator also may require many design trials in order to arrive at an 
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acceptable design for all the desired criteria. If there is a reason for 
using a nonlinear device as a compensator, the calculations are as easily 
accomplished as for the inherently nonlinear system. Furthermore » the 
shape of the nonlinearity necessary to give the required output response 
nay be determined. There are no general restrictions imposed on the de=- 
sign of nonlinear compensators . This is a distinct advantage over pre- 
sent methods of design of nonlinear compensators , The method may be 
applied to the solution of systems with multiple nonlinearities with the 
same ease , 

D) Continuous Data Control System. 

The design of continuous data control systems may be achieved by 
the use of the time domain matrix by the introduction of the approxima- 
tion of the fictitious sampler and hold. The approximation and the atten- 
dant error is discussed in section 4-6. Basically, the sampling rate 
must be sufficient and a fictitious hold introduced in order to produce 
a continuous input to the controlled system. A continuous-data compen= 
sation network may be selected by any of the standard design techniques 
and then introduced into the feedforward channel and its compensating 
effects evaluated by the time domain matrix method, A single loop con- 
tinuous system v;ith fictitious sampler and hold is shown in Figure 
The fictitious sampler and hold are usually inserted in the feedback 
loop in order to take advantage of the filtering action of the forward 
channel on the input signal. The use of lead and lag networks as com- 
pensating filters may be investigated directly in the time domain by this 
method , 

If it proves desirable to investigate the use of a compensator in 
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Elements 

Figure 5-V, Continuous Control System With Fictitious Sampling, 

the feedback channel, the introduction of a second fictitious sampler 
and hold will allow for compensation by a D matrix in the matrix equa= 
tions. Then the discrete compensator D found to be desirable may be 
approximated by a continuous network in the time domain. By techniques 
of approximation in the time domain a netv/ork can be synthesized to 

2 9 

yield a prescribed output for a prescribed input,'' The time necessary 
for the calculation of the digital compensator and evaluation of a suita=“ 
ble continuous network may be considerably less than that for a design 
carried out for the continuous data system by standard s-plane or fre- 
quency response techniques , 

E) Time-Varying Sampled Data Control Systems, 

The design of time varying sampled-data control systems is possible 
with the use of the method of time domain matrices. The standard me= 
thods of analysis and design for non-time -varying systems, that is the 
z-transform, and frequency response methods, are not applicable. The 
time varying system is represented by the G(n,k) matrix and allows the 
general design procedures of the preceding paragraphs to be applied. 
Furthermore, the design of a continuous data time-varying system is 
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possible through the introduction of a fictitious sampler and hold,, 

This approach offers great possibilities to the designer, in the analy- 
sis and design of systems of which the dynamics vary with some variable « 
An important example is the high performance jet aircraft, where the 
dynamics of the aircraft vary with altitude and therefore time<> A com- 
pensator may be introduced and a design evaluated. This compensator may 
be a continuous data or sampled data network , The limitations , advan- 
tages , and purposes of each type of compensator are essentially the saffe 
as those for the non-time varying systems discussed in the pieceding 
paragraphs . 

If the variation of the time varying element is relatively large, 
it often is impossible to compensate with a non-time varying compensator „ 
In this case it is necessary to use a time varying compensator such as 
a time varying amplifier, or perhaps a time varying network. The systems 
of this class are often termed adaptive systems. It is useful to consi- 
der systems with time varying dynamics (poles and zeros of the system 
transfer function), time varying gain, and time varying sampling rate. 
P-^search on these concepts has been carried out.^*^ The area of adaptive 
sampled-data systems shall be considered further in the next chapter. 

Finally, it is possible to use a time varying compensator in a non- 
time varying system. For example, it is possible to compensate a sys- 
tem for specified criteria b3'’ the use of a time varying amplifier. 

This possibility was illustrated by an example in the previous section. 
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CHAPTFR 6 



THE ANALYSIS AND DESIGN OF M)\?TWA COiNTROL SYSTE^fi 
6-1 . Introduction 

An important class of control systems are those for which para- 
meter values of the system change as functions of some independent varia=' 
ble during the period of operation. An important example is the change 
in the system dynamics of a supersonic aircraft with altitude. Another 
example is that of a chemical process where a parameter may change as 
a function of the ambient temperature. Usually, a fixed invariable com- 
pensator will be adequate over a restricted range of operating conditions. 
If the system is expected to operate outride this region, then some other 
form of compensation is necessary. The compensator often used, varies 
or adapts to the changed operating conditions. 

At the present time there is no widely accepted fundamental defini- 
tion of an adaptive system, although several are advanced in the litera- 

11 12 13 

ture.' ’ ’A system which changes or adapts a parameter of the con- 
trolled system to drive the actual performance towards the desired per- 
formance shall be considered adaptive . The two basic segments of an 
adaptive system are l) the identification of the dynamics of the con- 
trolled system, directly or by means of a related variable^ 2) the 
generation of an appropriate actuating signal for the controlled system. 
Identification of the dynamics may be accomplished , for example , by 
measuring the impulse response , by measuring the response to white noise , 
or measuring a related variable such as the output to a known input 
signal. The adaptive actuator may be a nonlinear, time-varying, or a 
digital device. For complex systems, it is common to use a digital 
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computer in which case the sj'’stera must be treated as a sainpled-data 
system, 

6-2, The Use of A MDdel in an Adaptive System, 

One philosophy of design of adaptive systems incorporates the use 
of a system model in the input channel as shown in figure 6-1 , The 
output of the model is then compared with the actual output and the 
error used to drive the system towards the desired output, iinother form 
of a closed loop system using a model is shown in figure 6-2, This con- 
figuration has been called a conditional feedback system, the feedback 
of a control signal being conditional on the result of the error between 
the desired and actual system output. The necessary conditions for the 
desired performance are determined by the model and the Adaptive compu- 
ter provides the actuating signal. The model may be a physical simula- 
tion or analog of the process or a mathematical abstraction manifested 
as a set of equations stored in a computer. The Adaptive computer is 
often a special purpose digital computer, but may be a nonlinear or time 
varying controller. 

The application of time domin matrix methods to the analysis and 

design of adaptive systems follows the methods discussed in the previous 

chapters. The investigation of a sampled-data conditional feedback 

system shall illustrate the, basic approach. Consider a basic system 

where =1, the adaptive computer is a direct connection, and the 

feedback H(s) = 2, as shown in figure 6-3. The sampling period will be 

nd 

one second and the plant is a type I, 2 order system for which the 
system matrix is set forth in Appendix A. The goal of the design shall 
be to obtain a system output response for shifting system dynamics which 
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Figure 6-1. The Use of A Model to Shape the Input Signal 




Figure 6-2, The Use of A ^bdel in A Conditional Feedback System 
I 
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Figure 6-3, A Sampled-Data Conditional Feedback System 

approaches the output response for the desired system. The desired res- 
ponse is present in the loop as the output of the model, X(s)„ This 
desired response is compared with the actual and the difference used 
to drive the output towards the desired response. For exan’ple, if 
a = 1 for the unvaried system, then the output for a step input would be 
c] = Col’jmn{o, .3679, 1.00, 1.40, 1.40, 1.15, .?94, ,798, « • .] 
where (c| is the transposed column vector. This response is shown on 
figure 6-4 as curve number 1, Furthermore, if the varying system mani- 
fests its variation in a shift in the pole a, over the range a s .50 to 
a = 2.0, the simple single loop system output would vary in the limit 
as shown in figure 6-4 as curve 2 and 3. For the conditional feedback 
system of figure 6-3, one obtains the following matrix equations^ 

e] = r] + X] - 2C] (6-1) 

c] = [g1 e1 (6-2) 

Then , one may write equation 6-1 as 
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1 > 3Co 
1.3679 - 2C, 

2.000 - 20 ^ 

eJ = ■” 20 , 

2.40 - 20 ^ 

2,15 - 2C^ 

1.B94 - 20^^ 

1.798 - 20^ 

For the system with the model in the limit at a = 2.0, one obtains. 
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This output response for the system with the model is shown as curve 
number 4 on figure 6-4 and the response at the limit for the pole a = ,50 
is shown as curve number 5, It can be seen, that the system with the 
model gives a response more closely approximating the desired than would 
the uncompensated system for either the limiting pole magnitude of & = ,50 
or a = 2,0. The simple conditional feedback system illustrated here, 
yields a definite advantage for varying parameter systems. If it was 
necessary to duplicate the desired response more perfectly, it would be 
necessary to use an adaptive controller in the feedback loop as sho\m in 
figure 6-2. In a sampled data system this would usually be a digital 
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network or D block. The design of digital networks shall be discussed 
in the next chapter. 

6-3 An Vdaptive Gain System, 

One possible method of adaption for a system is to set a parameter 

of the closed loop system dependent on a variable related to a perfor- 

mance index. One proposal advanced in the literature considers the varia- 

12 

tion of a compensator pole based on a measured figure of merit. Ttte 
calculation of the figure of merit, such as ITAS (integrated Product cf 
Time and Absolute Error), is usually accomplished by a special purpose 
computer and consideration of the signals as sampled-data usually follows. 
One form of an adaptive system would use the magnitude of the 
sampled error to control the gain of the system. There is no signal 
storage involved in a simple gain adjustment system and therefore the 
compensation improvement is limited. However, this system illustrates 
the possibilities of this approach. Consider the system shown in figure 
6-5. The input shall be considered a step function, and the amplifier 
gain a function of the error. Then, for the equations of the system 
one obtains" 




Figure 6-5 » Maptive System V/ith A Variable Gain 
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k] = 

[a] = 
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0 



(6-5) 

function of the sampled error 



р] = [a] e] 

с] = [o] p] 



(6-6) 

(6-7) 

For a system where K = 1 and the amplifier gain is directly proportional 
to the magnitude of the error v/ith a minimum gain of ,?0, one obtains" 



P] = 



and 



C] = 
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(6-8) 



(6-9) 



The response for the uncompensated system is curve 1 on figure 6-6, 
while the compensated response is curve 2, The compensated system has 
less overshoot, shorter settling time and also a somewhat greater rise 
time. For this type of system, the designer has the choice of the mini- 
mum gain and the function of error that controls the gain of the ampli^ 
fier. If the pole of the plant is expected to change within the range 
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a =1 to a = ,5, then the designer may evaluate the response at a = „5 
for the compensated system. The response for the compensated system is 
shovm on figure 6-6 as curve 3 and the uncompensated response is shown 
as curve 4. It can be seen that, again, the compensated system will 
have a smaller overshoot, less settling time and approximately the same 
rise time as the uncompensated system, 

llie change in the output response effected by this compensation 
scheme may not be satisfactory for many purposes. It is not the intent 
of this chapter to treat ejdiaustively adaptive systems, but rather to 
illustrate the application of the time domain matrix method to adaptive 
systems , Another povjerful approach would involve relating the gain of 
the amplifier to the derivative of the error. The derivative of the 
error can be generated from the error sample magnitudes by the method 
of backward differences. The use of the error derivatives to control 
the gain or any other system parameter will yield a more desirable res- 
ponse than that using the error magnitude , A system which uses the 
error derivative to control the sampling rate is discussed in the next 
section, 

6-4, An Adaptive Samrling Frequency System 

Sampled-data control systems usually have fixed sampling frequen- 
cies which must be set high enough to give satisfactory performance for 
all anticipated conditions. It is useful to reduce the sampling fre- 
quency v;henever possible in order to extend component life and allow 
time sharing of the digital components, particularly the digital compu- 
ters. It is usually desired to have an efficient sampler. That is 
over a given time interval, fewer samples are needed with the variable 
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Figure 6=6 o Compensated and. Oncoinpensated System Response- 



frequency system than v;ith a fixed frequency system while maintaining 

essentially the same response characteristics. A study of an adaptive 

system which varies the sampling frequency hy measui'lng a system para=” 

10 

meter has been accomplished.’ It was shown experimentally, that a samp- 
ler whose sampling period Is controlled by the absolute value of the 
first derivative of the error signal v;ill be a more efficient sampler 
than a fixed frequency sampler. 

.Analytical methods of investigation were not available and in order 
to compare the experimental results with calculated results, it was neces= 
sary to develop a method of investigating sampled systems with varying 
sampling rates « Fortunately, the time domain matrix method with a time 
varying system matrix my be extended to investigation of variable fre^ 
quency sampling. Consider the system shown in figure 6-7. The sampling 
frequency is controlled by a function of the error, and could be con- 
trolled by a function of the sampled error if this was the available 
error signal as in a radar system. 




Figure 6-7. Variable Sampling Rate System 
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As was shown in the previous section, the derivatives of the error may 
be obtained from the sampled error by baclcv/ard difference formulas, 'nie 
output response may be written as * 

C (k )] = |_G (n , k)] E (k)] (6-10) 

where j^G(n,k)j is the time varying system matrix 

and C(k)J and E(k)j are time varying col'Jimn matrices. That isj, with 
the period of sampling changing with time, one cannot correctly write 

cj = ,02,03,0^, • • • ^ (6~ll) 

where the output samples are equally spaced T seconds apart, ifeen the 
sampling rate changes, samples occur at various times and in this case 
k equals the number of seconds elapsed from the time origin. 

As an example, consider a system with two values of sampling period ^ 
T = 1 second and T = 2 seconds. Obviously, it is more efficient to 
have a sampling period of tv/o seconds when the error is changing slowly, 
therefore the first derivative of the error may be used to switch the 
sampling pe ’iod from one second to two seconds , Consider a step input 
and a system with a transfer function I 

Furthermore, for Illustration, the sampling rate shall be considered 
to switch at the first sample when the derivative of the error is zero , 
which is at the first overshoot peak. Considering line A of table 6-1 
which is the response for this system with T = 1 second for all time, 
one can expect the same values of response for the first three seconds. 
At the output peak overshoot the sampling period switches to T = 2 
seconds and one obtains for equation 6-10* 
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Ttierefores, using the values for the given in appendix A for T - 1 
second up to n = 3 and then for T = 2 seconds one obtains the folloAd.ng<, 



where the star indicates the T = 2 condition* 
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(6-14) 



Ttierefore, the output response is; 



!(k)] 



G(0) 

C(l) 

G{2) 

G(3) 

G(5) 

G(7) 

G(9) 



0 

.3679 

1.000 

1.400 

1.1465 

.7083 

.8931 



(6-15) 



It can be seen that there will be fewer samples necessary for this varia- 



ble frequency and it is therefore more efficient. The resnonse of the 
variable sampling frequency system is given in table 6-3 . and it is essen- 
tially the same as the fixed frequency system with respect to overshoot p 



time (seconds) 


0 


1 


2 


3 


5 


7 


9 


A Output-Fixed T 


0 


.368 


1.00 


1.40 


1.15 


.802 


.994 


B Output-VariableT 


0 


to 

o 


1.00 


1.40 


1.15 


.708 


.^93 



Table 6-1. Response V/ith Fixed and Variable Sampling Period. 
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rise time, and settling tLme . This subject is treated further in re- 
ference 10, where it was shovn experimentally that a reduction in the 
number of samples of tv/enty-five to fifty percent may be accomplished 



6-5. The Use of the Time-Domain Iktrix for Analysis and Design of Adap- 
tive Systems , 

Adaptive systems may be treated as time-varying control systems In 
almost all cases , Therefore , continuous or sampled-data adaptive sys- 
tems may be profitably investigated by the use of time domain matrix 
methods. In fact, there may be no other method than can be used except 
solution by numerical methods ultilizing the digital computer. 

One of the necessary steps in the adaptive process is the identi- 
fication of the system, that is, the impulse response or the transfer 
function, With the use of time domain matrices it is possible to de- 
termine the pulse response when the system is at rest as outlined in 

chapter two. If the system is not at rest at any time, then the iden- 

13 

tification of the dynamics must be accomplished by some other means, 

A predictor or learning system may be designed by the use of back- 
ward and forward differences , Tlie evaluation of a forward difference 
equation allows the system to predict its next few values and adjust 
for optimum conditions , For example , the Gregory-Newton forward inter- 
polation formula may be programmed in a digital network for the pre- 
dicted output as I 



'n+1 



3c - 3c , + c ,, 

-^n n-1 n-2 



i-V os ^ 

(6-jLo / 



using the present value and two past values requiring two storage ele-^ 
ments , The accuracy of prediction increases with an increase in the 
number of storage elements and if three past values were to be used 



in the calculation, then the formula to be programmed would be I 



c T =4c “6c ,+4c --c , (6-17) 

n+1 n n-1 n-2 n-3 

On the basis of the predicted value of the output response, the adap- 
tive system may adjust a system parameter or add a signal to drive the 
actual output towards the desired output. The accuracy of this method 
improves with an increasing sampling frequency as would be expected. 

This adaptive system using equation 6-16 and 6-17 may be implemented 

by the use of digital logic networks or operational amplifiers and el- 

P 

ectronic samplers , Equations 6-16 and 6-17 may be written in matrix 
form to aid in the calculation of the predicted values . 

The use of a special or general purpose computer in an adaptive 
system can be investigated in general by the method of time domain 
matrices , One adaptive system uses the integrated Product of Time 
and Absolute Error as the figure of merit and obtains the magnitude 
of this figure of merit by calculations on the sampled signals. These 
calculations are accomplished by a programmed numerical method and can 
be written in matrix form. Furthermore, the method of steepest des- 
cent may also be calculated by the matrix approach. In general, the 
use of numerical methods in the calculation of the figure of merit and 
adjustment necessary in adaptive systems may readily be accomplished 
through the use of time domain matrix methods , 
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CHIPTER 7 



THE DESIGN 0? DISCRETE-DIT\ COMPENS ITORS 
7-1. Introduction 

In the design of control systems it is desired to satisfy a set 
of specifications and response requirements. In many cases, this re- 
quires the introduction of a compensator in the control loop. The de- 
sign process may be carried out by the use of the time domain matrix 
for sampled-data systems and for continuous data systems by the intro- 
duction of the fictitious sampler. 

The introduction of a discrete-data compensator in the control 
loop will allow the designer to achieve a required output response. 

The com)-)ensator operates on the discrete-data input and yields a dis- 
crete-data output, Tfiis compensator may be realized by a special or 
general purpose digital computer, a logic network, or a circuit com- 
posed of operational amplifiers and electronic samplers. 

The process to be controlled is usually a continuous process and 
therefore it is necessary to investigate the response of the system 
between sampling instants while carrying out the design. Design ly 
z plane or root locus methods does not include this possibility and 
often results in unsatisfactory intersample system response , The use 
of time domain design techniques has two basic advantages I 

1) the design is carried out directly in the time domain 

2) the intersample responvse is accounted for in the design, 

A sampled data system with a discrete data compensator is shown 
in figua-e 7-1, The D (s) is the sampled data transfer function of 
the compensator. 




Figure 7-1, A Sampled-Data Systers with A Discrete-Data Compensator-. 



7-?. The Physical F.ealizability of Discrete-Data Compensators, 

Trie digital compensator operates on the input signal and yields 
a transformed discrete tine sequence output. The z-transform equation 
for this operation may be written for figure 7-1 as * 

P(z) = D(z)E(z) (7-1) 

Then, it can be seen that this operation may be wi’itten as a matrix 
equation as follows I 

P] = [d] e] (7-2) 



where 



[dJ. 



^0 

d. 



0 0 
do 0 
d, do 



0 

0 



n 



when the compensator is time invariant. If the digital compensator 
is time varying, then one may write* 



p] = [D(n,k)J e] 



where 








t 



d(0,0) 

d(l,0) 

d(2,C) 



0 



0 

0 



[D(n,k)] = 



d(l,l) 
d(2,l) d(2,2). 



It is necessary to determine what restrictions are imposed upon the D 
matrix by the requirement of physical realizability, 

A system is said to be physically i-ealizable , if the output signal 
of the system does not depend upon future information of the input sig- 
nal . In the matrix equation 7-1, this is expressed as the requirement 
that all the elements of the D matrix above the main diagonal be zero. 
The truth of this statement can be seen from the follo;<dng matrix equa- 
tion. 



P] = 



do d^ 0 



di 



a 



0 

0 



d, 



a 



e] 



(7-3) 



Expanding the second row (n = l), one obtains ^ 

p, = d,eo + doe, + d^^e^ (7-4) 

Therefore, p, would depend upon the future input signal which is not 
available at the second sampling instant (n = l). Therefoi'e, it is 
necessary for d =0, and all elements above the main diagonal to be 
equal to zero. 

Furthermore, for a stable compensator, it is necessary for° 

d 



limit 
n -»• 00 



n+1 



d 



n 



N< 1 



(7-5) 



If this limit is equal to one then the D block is a steady-state oscil- 



F9 



lator. 



4^..i 




The discrete data compensators may be realized byt^ 

1 ) Digital Programming 

2 ) Delay Line Ketwoi-ks 

3 ) Discrete-data RC networks 

4 ) Analog Computer Elements . 

7-3. The Sensitivity Piatrix of the Sampled Data System with A Discrete 
Compensator . 

A logical quantative measure of the control property of a feedback 
system is the sensitivity which is defined as the relative change in 
the system transfer function T, divided by the relative change in the 
plant In the z transform notation, this my be expressed as* 



0 - oT/ X _ dTCzVlCz) 

^ “ ^G^^^ ■ 'dG{zS/G(zT^ 


(7-6) 


where T(z) = “^{ 7 }- 




For figure 7-1, one obtains! 






(7-7) 


and ^ 

^ = 1 . GTiTDfzT 


(7-8) 



This system possesses two degrees of freedom, which permits indepen- 
dent realization of T and S by means of the G and D components. Equa- 
tions 7-7 and 7-8 may be written in matrix form as * 

[t] = [g] [d]{[i1*[g][d]|-1 (7-9) 

""" [s] -{[i] . 

Then the sensitivitj’’ specifications fix the sensitivity matrix and 
therefore , 
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I ' •'i 

p ■ 



UP 









Ui 



(7-10) 



[o][d] = [S]-^ -[I] 

The required output specifications fix the T matrix and since J 

[t] =[g] [d] [s] (7-11) 

one obtains! 

[g] =[t] [s]-1[d] (7-12) 

Therefore, equations 7-12 and 7-10 may be used to obtain the required 
system matrices G and D. If the plant G is fixed by power and load 
considerations, then a new compensator, D, , must be introduced in the 
control loopj. 

The sensitivity of an uncompensated and the compensated system is 
sho’/m in figure 7-2. The uncompensated system is a type one, second 
order system with a zero order hold. The compensation reduces the over- 
shoot for a step input, to half the uncompensated value while reducing 
the rise time by half a sampling period. The magnitude of the sensiti=^ 
vity response of figure 7-2 is greatest during the first three sampling 
periods and settles out most rapidly for the compensated system, 

TTie use of the sensitivity matrix or response, aids the designer 
in understanding the effects of compensation by a discrete data network. 
Furthermore, the variation of the system response with a system changing 
with time may be studied for adaptive systems by considering the sensi” 
tivity. 

7-4, Design of a Closed Loop Discrete Compensator by Means of An Open 
Loop Discrete Compensator. 

For the single closed loop samrled-data system shown in figure 
7-3, a design of a discrete compensator is usually required in order 
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Figure 7-2, 1116 Sensitivity of the Compensated and Uncompensated 

' System 
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7-3. Error Sampled System with A Discrete Compensator 

for the system to meet the necessary specifications. The output res- 
ponse at the sampling instants may be written in a time domain matrix 
equation as I 

This complex relation involving D causes the investigator to consider 

2 

the design of an open loop compensator as shovm in figure 7-4. 



R(s) 




7-4. An Open Loop Discrete Compensator 

Then, the matrix equation for the sampled output response my be 
written as ’ 

C] = [g] [Do] r] (7-14) 

If an open loop compensator may be determined , then a closed loop com- 
pensator may be found , To find the relation between the open and closed 
loop compensator, equations 7-13 and 7-14 are equated and solved for D 
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yielding*. 



[d] = W{[l] -IDo][o]|-1 



(7-15) 



It can be seen that this method takes the plant G into consideration in 
the determination of the closed loop compensator. 

An example will illustrate the design procedure . Consider a sys- 
tem with a transfer function G(s) = ^ zero order hold, 

and a sampling period of one second , A design of a compensator shall 
be accomplished for the best compromise for a step and ramp input. The 
open loop compensator is chosen to yield a system output with a final 
value of one for a step input. Then, one has' 



- 1 > 



Do(z) = (1 + d,z"-") (1 



) ( 



1 + d, 

-1 



) 



= ( xVd^ - d,z“^) 



(7-16) 



This compensator will have a minimum complexity since the use of further 

-2 

delays such as d^z ', increases the complexity of the closed loop com- 
pensator. The final value in the time domain, for a step input, yields 
a value of unity for the system output as expected . The final value 
theorem is written and evaluated as follows 



limit 
t CO 


C] 


lim 
t-* oc 


[g] [d.] r] 






(7- 






_ 














0 


0 0 




1 


0 


0 • • . 


1 




gi 


0 




(d, -1 ) 


1 


0 


1 


limit 


c 


• 


^ 1+d, ^ 


-d, 


(d,-l) 


-1 


1 


t-* CO 


« 

0 


• 

• 


0 


-d, 


(d, -1 ) • • 


1 






^n+1 


1 


• 

1 

• 




• 

• 










i ; 






• 


- 



limit 
n — cc 



in±Llf:in 

1 + d. 



= 1 
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The output response for the sampling instants may be written for a step 
input as I 



=] = <1^7’ 



0 0 
gi 0 



0 

0 



g2 gl 0 
g3 ^2 





10 0 0 


1 




(d,-l) 10 0 


1 




-d, (d,-l) 1 0 


1 




0 -d, (d,-l) 1 


1 

. 




r 

o o 

1 


• 






0 0 
gi 0 



0 

0 



g2 gi 0 
g3 g2 gi • • 



1 

di 

0 

0 



(7-18 ) 



Therefore, the output at the sampling instants is ° 

, n >, 1 



°n = 



e + d, g , 
^n ’ ^n-1 



(7-19) 



1 + d, 

The response between the sampling instants may be evaluated from° 

C(m)] = [c(m)] [Do] r] (7-20) 

Then, for half way between the sampling instants (m = 1/2), one obtains 
for the output response" 

S. 



c (1/2) = fn-1/2 ^n-1.^, 



n 



For a ramp input, the output response may be written as' 

where r] = column ^0,1, 2, 3>4, • • ‘I 
Then, one obtains t 



(7-21) 



(7-22) 



Co = c, =0 
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c 

n 



^n-3 

1 + d. 






=n-3 



n 2 (7-23) 



The output response is then easily evaluated for the step and ramp in- 
puts and the results for this example are shov/n on figures 7-5 and 7-6. 
The designer would be able to choose a value of d on the basis of the 
specifications for a step and ramp input. In this case, as a compromise, 
the designer might choose d, = -.60. Then, the matrix for the closed 
loop compensator D may be evaluated from equation 7-3,5. Since the sys- 
tem G(s) is second order, the closed loop compensator required will be 
third order. 

If more control over the output response is desired by the designer, 
then one may use a higher order Dq . For example , a higher order compen- 
sator might be * 

With this form of compensator, a deadbeat response to a step input is 
possible. In general, this method of compensation in the open loop, 
gives the designer strong control over the system output response . The 
method may be applied to linear sampled-data and continious systems, 
and time varying sampled-data and continuous systems. 



7-5. Design of The Closed Loop Compensator By Time Domain Evaluation. 

A simp3e method of the design of the closed loop discrete compen- 
sator should not be overloolcedj that is, the arbitrary choice of a D(z) 
and an evaluation of the output response by the use of time domain ma- 
trices . The evaluation of the output response between the sampling 
instants avoids the possibility of hidden oscil3.ations . The choice of 
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Figure 7-5 • Response of the System to A Step l^put 




Figure 7-6, Respoase of the System to A Ramp Input 



the D(z) may be aided by all the standard z-transforra techniques ^ and 
evaluated by the time domain matrix. For the example, of the previous 
section, one might choose on the basis of z -plane techniques a conpen” 
sator as * 



D(z) 



K(z - .367?^?) 
(z + d) 



(7“25) 



Then, evaluating the output response for K = 2, d ■= ,50 and d = ,72, 
as shown in figure 7-7, the designer would choose the compensator on 
the basis of the specifications. It is interesting to note that for 
K= 2, d = .72, one obtains 



G(z) D(z) = 



(. 3679) 
(z - l) 





z + 




- » 7.358 

- IrroT 



(7-26) 



This type of compensation in the z plane, conuionly called cancellation 
compensation is misleading since the z-plane only accounts for the sampl= 
ing instants , The total res’^onse for all time may be evaluated using 
the time domain matrix. If Kis set equal to 2.718 and one obtains" 



G(z) D(z) = 



(7-27) 



and the closed loop, z plane response is* 

C(z) = z“^ R(z) (7-28) 

The response for this system is shown on figure 7-7 as curve number 3. 

It is obvious that it is necessary to consider the intersample response 
in most design problems . 

7-6. The General Characteristics of the Closed Loop Discrete Compensator, 
It is worthwhile to look at the general form of the discrete compen- 
sator D, The compensator may be written in the z-transform ast 
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(7-29) 





The time domain matrix is written as ° 



1 0 0 0 0 



0 o o 



d, 1 0 0 0 



[d] = K dj d, 1 0 0 

d3 dj d, 1 0 



(7-30) 



For a first order compensator one has ’ 



D(z) = 




(7-31) 



where the zero is in the right half of the z plane , almost always the 
case , Then , one obtains ° 



Therefore, if the values of d, and dj are obtained as necessary for com- 
pensation, the necessary D(z) may be realized as equation 7-31, If it 
is necessary to specify four elements of the D matrix then the compensa- 
tor must be second order. That is, if d, , dj d 3 » d^ are specified, then 
D(z) must be of the form' 



The sampled impulse response of the discrete network weights the 
past and present values of the input signal, yielding the output signal 




(7-32) 



For this compensator, the matrix values are; 

d, = -(a+b) , dj = b(a+b) , d 3 = -b (a+b) , 



(7-33) 




(7-34) 
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at the sampling instants , 'There are three possible sampled inpulse 
response time series for a stable discrete compensator. These possi- 
bilities are illustrated by figure 7-? which also lists the necessary 
first order compensator and its z-plane pole arid zero locations. An 
unstable discrete compensator occurs when b > 1 in case c, where then 
the alternating series is a diverging time series . 

It is worthwhile to investigate further, the alternating weighting 
series of the compensator of case c, Tliis discrete compensator would 
normally be used to stabilize the system while maintaining a desired 
rise and settling time . From equation 7-32 one can observe that if 
(a+b) is greater in magnitude than one, the second term of the con- 
verging series is greater than the constant term of one. Therefore, it 
is actually possible to v/eigh past sampled inf ormti on with a greater 
value than present information. This condition my be used, although 
usually to weigh past data with the greater value would have an unsta- 
bilizing influence . A picture of the sampled impulse response of the 
discrete compensator is useful in designing a compensator as shown in 
the next section. 

7-7. Time Domain Design of A Discrete Compensator. 

Hie design methods discussed in the previous sections depend either 
on standard techniques on the z-plane or on the choice of a trial com- 
pensator and evaluation in the time domain. These methods are powerful 
and are the ones used predominantly in practice. However, it would be 
of value to the designer if one could design the compensator directly 
in the time domin, Tvro methods of compensation design in the time do- 
main by the use of time domain mtrices have been developed by the author 
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a) D(z) = 9 l>|b[>|a[ 




Figure 7-8 * The Stable Sampled Impulse Response of The 

Discrete Compensator 
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and are presented in this section 



A. Synthesis at the Sampling Instants. 

For the error compensated single loop sampled-data system shown in 
figure 7-3 » one my write the follov/ing time domin matrix equations' 



C] = [g] p] (7-35) 

P] = [d] e] (7-36) 

’^here ^ 

Also, the equation for the response between the sampling instants may 
be written as I 

C(ra)] = [G(m)] p] (7-38) 

Tlierefore, if the designer specifies the values of the C mtrix in equa- 



tion 7-35 on the basis of the required response, the values of the P 
matrix nay be determined. Then, with a known input R, and a specified 
C matrix, the error matrix may be evaluated from equation 7-37, Since, 

P and E are then known, the required D matrix may be evaluated from 
equation 7-35. It must be noted that the designer has no control over 
the intersample response when using this method , When the P matrix is 
found on the basis of the desired output C at the sampling instants, 
the output betv/een the sanmling instants C(m) is then calculated from 
equation 7-38 . 

Once the values of the sampled output are specified , the P matrix 
may be determined from equation 7-35, or by vn?iting* 

p] =[g]"^ c] (7-39) 

It is usually convenient to use equation 7-35 and avoid the inversion 
of the G matrix. It is possible to specify an output response at the 
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sampling instants which the compensated system will be unable to yield 9 
in which case the calculated compensator, D(z), will be unrealizable. 

For most design problems, it is sufficient to find the first few ele- 
ments of the D matrix and determine a D(z) compensator on this basis 
from equation 7-32 or 7-34. 

An example will illustrate the procedure and the possibilities of 
conj^ensation. Consider the system of figure 7-9 which has been considered 
in the previous sections. The output response of the uncompensated sys- 
tem is shown as curve number 1 in figure 7-16. 




Figure 7-9. Compensated Sampled Data Control System. 

It is desired to lower the rise time (for 90 ®/o of the final value) 
while lowering the maximum overshoot to 10 percent. Then, on this 
basis, set I 

c] = column ^ 0 , .50, 1.0, 1.10, • • (7-40) 

and using equation 7-35 evaluate the necessary values in the P matrix. 
For the first four sampling instants, one may v/rite ° 
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(7-41) 



Co 




~ 0 


0 


0 


0 


Po 


Cl 




.3679 


0 


0 


0 


Pi 


C 2 




.7675 


.3679 


0 


0 


P 2 


C3 _ 




.9145 


.7675 


.3679 


0 


P3 



Using this matrix, the may be evaluated step by step as follows' 



Po 

Pi 



= ■ ■ 1 = 1 '3KQ 

.3679 .3679 

_ 1.0 -(.7675) fl.359) 
g, .3679 



-.1170 



(7-42) 



P2 



“S 3 P 0 “S 2 P 1 



-.1441 



Then, the D matrix my be evaluated on the basis of equation 7-36 and 
the calculated error values as follows ' 



p] = [dJe] 



1.359 




do 


0 


0 


1.0 


-.1170 


= K 


di 


do 


0 


.50 


-.14a 




d2 


di 


do 


1 

0 

0 

• 



Therefore, the d values arej 



(7-43) 



Kdo = 1.359 , Kd, = -.7965 , Klj = + .2542 
and the D matrix may be written as ; 



k[d] = 1.359 



10 0 

-.5^61 1 0 

+.1870 -.5?61 1 



(7-44) 



If a first order compensator is used to realize this compensator, equa- 
tion 7-31 through 7-33 will apply. Using equation 7-33, the compensator 
is found fromt 

(art) = .5861 , = .3191 

Therefore, a = .2670 and the compensator D(z) my be written as' 
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(7-45) 



= 1.359 (1 -.5S61z”^ + .1870z 



.0597z"^ + ,0191z 



For this compensator, the total response nay be evaluated using equations 
7-35 and 7-38. The response of the uncompensated and compensated sys- 
tems are shovm on figure 7-10, as curves number one and two respectively, 
The output response between sampling instants is not directly taken in- 
to account in this design approach. Therefore, if too stringent require- 
ments are imposed on the response at the sampling instants, the inter- 
sample response may become unacceptable. The designer will usually eval- 
uate the intersample response for a design choice by using equation 
7-38. For example, if the specifications called for a maximum overshooxi 
of five percent for the previous example, one might attempt to set C3 
equal to unity rather than 1.10 as was previously done. This design 
change will not effect the values of po and p, or do and d, that have 
been determined . Therefore J 

Po = 1.359 , p, = -.1170 

C3 "g3Po “gaPi 

and P2 = ^ = -.4159 (7-46) 

Then, the gain and pole and zero of the first order compensator may be 
determined as follows I 



= Pa “.5K3, = -.4159 -.5(-.7965) = -.01765 



(7-47) 



where K = 1.359 and Kl, = -.7965 as previously. 



Therefore, from equation 7-33, 



(a+b)=|d,| =,5861 



and 
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u 



d 



b = 






-.01299 

75861 



= -.02216 



(7-48) 



a = .6083 

Then , the compensator may be vrritten as ° 



(7-49) 



= 1.359 (1 -.586lz"^ -,01299z“^ -.000287z"^ - . . .) 

The output halfway between the sampling instants may be calculated from 
equation 7-38 as I 

(7-50) 



7i/2)] = 



[g(1/2)] 


[“] E, 


] 




.1065 


1.359 




.6166 


-.7965 


e] 


.859 


-.0177 


.9481 


-.0004 




.1447 


1 




.145 


.7531 


.5 


z: 


.8255 


.6744 


0 




1.0509 


.5933 


0 




.9305 



This response would have an overshoot of five percent between the sample 
ing instants and the design would marginally meet the specifications. 
This response is shown as curve number 3 on figure 7-10, If the desig- 
ner chose to set C3 equal to 1,050, then the complete output response 
for a compensator calculated for this new set of response values at the 
sampling instants is shown as curve number 4 in figure 7-10, The maxi- 
mum overshoot for this compensation is 6.5 percent. A comparison of 



the rise time, maximum overshoot, and settling time for the uncompensated 
and compensated systems is presented in table 7-1 . The designer would 
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hi 

u 



•:3P .ond U)\;^ ".’■ap- :;;sa'*- E®?por 







0 ’ 






I 




chose the compensation that ’^rovided the closest approximation to the 
original specifications. It can be seen that it is a povrerful method 
of design of discrete compensators , The compensator may be located in 



Curve Number 


1 


2 


3 


4 




Mixiraum Overshoot 


45 


10.0 


5.0 


6.5 


90 Rise Time 


1,S2 


1,65 


1,65 


1,65 


seconds 


2,5 % Settling Time 


16,0 


4.S 


6.4 


5,2 


seconds 



Table 7-1 , Indices of Response For the Compensated and Uncom=» 
pensated Systems, 

the feedback channel or in one loop in a multiloop system and the design 
method loses none of its usefulness. However, if the designer has a 
set of stringent specifications , it may be necessary to account for^ the 
intersample response in the design. 

B, Synthesis of the Discrete Compensator Accounting for the Intersanrple 
Response . 

It is possible to account for the response between the sampling 
instants hy the use of the time domain matrix equation for the inter- 
sample response which may be wj’itten as I 

C(m)] = [c(m)][D] e] (7«51) 

Then, the designer may chose the first few elements of the D matrix by 
considering the sampled output v;hile simultaneously considering the re- 
sultant intersample response value . Therefore , for example , the desig- 
ner will chose the Kdo value of the D matrix by examining the resultant 
c, and c ^ values at the same time. This process is carried on step by 
step until all necessary values of the compensator matrix are determined , 
In this manner, the next step would invola’e the determination of 

by a choice of c, and Ct c. The choice of c and where is 

< 1 , p n n-A 



no 



m 




ucually equal to one-half to yield the response nidway between samples ^ 



is facilitated by the determination of an expression relating the two 
variables and its graphical presentation. The three basic matrix equa= 
tions are° 

c] = [g][d]eJ=[d][g]e] (7-52) 

C(m)] = [G(m)J [d] e] = [d] [G(ra)] e] (7-53> 

e] 3 r] -c] (7-54) 



where order of matrix multiplication may be interchanged for time in- 
variant systems. Then, for the typical system with g© = 0» and a unit 
step input , one obtains “ 

c, = Kdogi (7-5:') 

c^5 = K3og^5 (7-56) 



It is usually necessary to choose Kdo such that the value of the first 
sampled response is less than the final value in order to achieve small 
maximum overshoot . In the previous example , for instance , Kdo was 
set equal to 1,359. Then, one my solve for c ^ in terms of c^ , as^ 

(7-57) 




Equations 7-55 through 7-57 allow the designer to choose the value of 
Kd(j while taking into account the intersaraple response , When the valce 
of Kdo is chosen, the value of e, is determined and then since 
d (3 = 1 by the definition of equation 7=79, one obtains" 





K 


0 


1 






gi 






K3,, 


g2 





(7-53) 
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Cii^i^ P*' 






and 



(1/2)J 



K 


S.5 


1 


Ki, 


^1.5 





(7-59) 



Therefore , 

C 2 = Kd,g, + K(g2 + g^e, ) 

and 

°1.5 = K(gj.5-*e.3*.) 



(7-60) 



(7-oi/ 



where the g <, K, and e< are known. Solving for c, in terms of e, 

n l.v * 

one obtains the general form of ° 



'^ 1.5 



(— ) 

' gi / 



c, + Q 



(7-62) 



where Q is a constant. In general, the linear relation may be written 
as ° 



'n-A 



( 



^1= 

"iT 



■) c„ + Q„ 

n n 



(7=t.3) 



It is often useful to show this linear relationship in the form of a 
graph. 



Consider the example of the previous section A, where the specified 
maximum overshoot was five percent. The gain constant K shall again 
be set equal to 1.359, so that c, = .50 as in the previous example. 

The design procedure in section A did not account for the intersample 
response and it was found that this resulted in the responses presented 
in figure 7-10. The design of section A to satisfy the overshoot re- 
quirement would have a response as in curve number 4 on figure 7-10. 

Now, in order to carry out the design while accounting for the inter- 
sample response, one obtains the equations for Cj and c, g terras of 



112 



^ f ^ 








m .^- 

( ^ 

\ 




Kd^ from” 



and 





1.359 


0 


1 


c] = 


Kd, 


.3679 


.50 






.7675 


62 



(1/3)] = 





1.359 


.1065 


1 


= 


Kl, 


.6166 


.5 




Kd2 


.859 


ej 



Therefore . one has i 



C2 = .3679Kdi + 1.2931 



and 



'1 ^ 



.1065Kdi + .9104 



Then ^ is related to Cj as! 



(' 7 -^ 4 ) 



(7-65) 



(7-66) 

(7-67) 



5 = .2?96c2 + .5360 (7-68) 

and this relation is shovm on figure 7“11. Considering the equations 

7-66 through 7-68 and the figure 7-11, the designer might choose 

C2 = .980 and c. = .8234 as a cocnromise. Then, one finds Kd. = -.851, 
-I- .5 

and d^ - -.626. The next step is to wi’ite the equation for C3 and 
^ noting that the maximum overshoot for this design vdll occur in 
this interval. One obtains from equations 7-64 and 7-65 J 



3 = .3679Kd2 + .9648 


(7-69) 


• 

^ ^ = .1065Kd2 + 1.0194 


(7-70) 


2 ^ = .?896c3 + .740 


(7-71) 



Equation 7-71 is plotted on figure 7-11, and this curve shows that the 
overshoot may be limited to five percent. If C3 is chosen as 1,040;. 
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m 



then = I 0 C 4 I and Kdj = ■>•.2045. A compensator is then determined 
from equation 7-33 as ' 




(7-72) 



The total output response is then determined and is plotted in figure 
7-12 as curve 5. The uncompensated system response is again curve 
number 1 , and curve 3 is the response for the design of section A 
which did not account for the intersample response. The rise time and 



been reduced 50 percent and there is no undershoot present. 

As another example, consider the single loop system as shown in 
figure 7-9 i. where the sampling period is again one second, the input 
is a unit step, and the system has a third order transfer function. 
The transfer function and zevo order hold are written as I 



and then the system transfer matrix is given in appendix A, table A-1 
section II, The uncompensated response of this system is unstable, 
and the designer might try reducirig the gain of the system in order to 
achieve stability. If the gain is reduced to two, or one-half of the 
unstable gain, the output response is shown in figure 7-13. This system 
has a mximum overshoot of 75 percent and a settling time of greater 
than 10 seconds. If the designer desired to maintain the same rise 
time while reducing the overshoot to 20 percent, it would be necessary 
to introduce a compensator. Using the discrete compensator design 
procedure, for a system gain of 4 as originally specified, one writes' 



overshoot are equal for each design, while the 2,5 % settling time has 




(7-73) 



s‘^(s-t-l) (s+2) 



114 




Figure 7-11 , The Relation of the Litersample Response and the 

Sampled Response 
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01^01 




0 ] 





K 


0 


1.0 




Kd, 


.3362 


ei 






1.1869 

1,6736 





and 



( 1 / 2 )] 





K 


.0582 


1.0 






.7845 


ei 




Kda 


1 ,4789 


_ 



! 7 -'’ 4 ; 






Choosing K = 1,0 and therefore using the original system of 4» one 
obtains ° 

c. = .3362 



and 



,6638 



(7=7o; 



Then, one uses equations 7-74 and 7-75 to obtain the equations for 
and the relation between them, as° 



‘'2.0 


= ,3362d, 


+ 1.410 




^ 1.5 


= .0582d, 


* .8231 


07-7^) 


*'1.5 


= .17320 2 


+ .5789 


^7-^9; 



Equation 7-79 is plotted on figure 7-14, Now, choosing approxiinatex;y 
equal to 1,0, one might use d, = -1,0, Tnen, obtaining the relations 



and - one has^ 






2.5 






c, p, = .3362dj + 


1.0266 


t7.80; 


C- ^ = ,0582d2 + 

2,5 


1.172 


(7-81 j 


c, ^ ^ ,1732c^ + 

^ 0 y 


.9942 


1,7-32; 



Equation 7-82 is plotted on figui'e 7-14 and from this figure one may 
choose Co approximately equal to ^ about 17 percent oversnoot. 
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A first order compensator 



Then it is found that may be set at ,5, 



is determined for these values and one obtains' 



D(z) = 





The response for the compensated system is shown in figure 7=»13o Tne 
maximum overshoot is 20 percent » the rise time is reduced to I*? seccrds 
and the settling time is reduced to 6,4 seconds, 

Tnese design methods are equally applicable for other system input 
signals. This will be shovm by designing a compensator for the system 
considered in this section subjected to a unit ramp signal Iripat, Again 
the sampling period is one second and the system and hold transfer fua. 
tion is ° 



The input signal and the uncompensated system response is shown on 
figure 7-14, The first two samples have a zero values, ®ay 
write for the output response ° 



Then^ choosing K = 2 so that Cj = .735? ^ one may write for tne next 
sample interval' 




s^(s + 1) 



(,7-8a) 



0 0 




K .3679 1.0 

M, .7675 e, 

.9145 63 



and 




K .1065 0 

Kd, ,6166 1 .0 

Ms .859 ej 



(.7=80) 
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3679M^ ^ 2,465 



(7-87) 



o 



5 = ,1065K3, -»• 1.5026 



(7-88) 



and 



2.5 



2896c3 5- .7977 



(7-89) 



From tnese equations and the nature of the responses one might set 
Kit ~ ”o50. Then repeating these steps fcr the next interval, one 
might let = * ,20. The value of Kd^ chcsen would depend largely 
on the overshoot limitation. In this case, tnere would be no overshoot 
present. The total response is showi in figure 7-14 as curve number 
two. The steady state error of the compensated system is sixty percent 
of the steady state error of the unccnpensatcd system. The compensator 
to yield this response is found to be! 



If it is required to have a shorter rise time, while permitting an 
overshoot less than five percent, it is necessary to raise the value 
of K to three. Then if the usual steps axe carried out, one obtains! 



Then, the total response is found to be as C!.irve number three in figure 
7-15. The rise time is greatly reduced, while the maxim’um overshoot is 
less than two percent, :\lso, the steady state error is reduced to 42 
percent of that for the uncompensated system. The compensator neces- 
sary to realize this response is found to be! 



If the input signal was expected to alternate between a ramp and step 
signal, it would be necessary to design a compromise compensator with 




(7-90) 
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Figure 7“15» ' System Response to A Ramp Input 



1 




respect to the design requirei«sntt „ 7ni> is necessary, since a con- 
pensator r.c be cptinun '^ith respecr “ ; c set of conditicns, can- 

t 

not be expected to be crtiJEcai for a cr.ingec. set of operatirjg conditior»s 
and requirements o 

It can be seen that the design opr-rcacbes presented in this sec- 
tion allows the designer a wide latitude cf the choice of specifica- 
tions, and affords direct control over toe output response,, Furthermor 
the time domain approach gives the designer a direct picture of the 
time response, therefore avoiding Insxoura^e and unwieldy correlation 
theorems for the z plane „ Fne desig?. method which accounts for tne 
intersample response allows the designer r.c wield control over the tots- 
time response. It is obvious that this design procedure nay be pro- 
grammed in a digital computer, the designer supplying the specifica- 
tions and system transfer matrix as the i^^put to the computer, 

A discrete compensator may be determiac-d for a continuous system 

which is approximated by a fictitcue sa.Tfpj.er and hold as discussed 

previously. Then the calculated D(s) may be sy-nthesized bj*- a continue 

9 

ous compensator in the time domain. Therefore, this technique is 
not limited to sampied-data systems , but ira,y be applied to any type 
of continuous data system. 

The location of the discrete ccm-penaatcr is not limited to the 
error channel and my be located In <any arbitrary loop in the system 
as feedback compensation. 

Therefore, two methods cf design, directly In the time domain, 



have been presented. These design procedures are entirely flexible 
in appli jaticn, accurate in calculation, and rapid In solution. The 



design in the time domain gives the designer a complete insight into 
the total time response of the system^ No other existing design pro- 
cedure can be carried out directly in the tic» domain, nor can any 
existing design procedure provide the flexibility and accuracy of the 
time domain matrix method. 
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CHAPTER 8 



COI« FUSIONS 

P-1. SuMnary of Results 

The aim of this dissertation was to present a new method of 
engineering analysis and design for c^anrplex control systems. This 
method is the time domain infinite aaatrix method , The formulation 
of the infinite matrix follows from the iSonvolution summation of 
sampled data systems. The mathematical basis of the time domain 
matrix formulation is presented in a discussion of the applicable 
concepts of infinite matrices and sequence spaces , This method of 
analysis and design is applicable to both continuous data and sampled 
data systems. For continuous systems it is necessary to introduce 
a fictitious sampler and hold of sufficient sampling rate to effect 
an accurate approximation. 

It is possible to analyze and design linear, nonlinear, and 
time varying systensof the continuous or sampled data class. Sam- 
pled data, time varying systems may not be investigated by any other 
existing method. Furthermore the investigation of nonlinear systems 
is greatly simplified by the time domin approach, fiiiltiloop sys- 
tems may be treated with ease and the signals at intermediate points 
throughout the loops are readily available. Also, systems with 
multiple nonlinearities may be investigated, for which there is not 
a presently available method of analysis and design. 

Two methods of design of a discrete compensator for a sampled 
data system are presented , These methods are accomplished directly 
in the time domain and allow for a compromise of specifications in 
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the time domain,, Also the response between sampling instants is 
accounted for in one of the two design procedures , 

The time domain matrix method iaa .7 be readily programmed on a 
digital computer and therefore provides a rapid analysis and design 
technique . 

8-2. Further Conclusions 

It is an important advantage that the design and analysis of 
systems using the time domain matrix method takes place directly in 
the time domain and not in any transformed complex variable domain. 
This advantage aids the designer in understanding and controlling 
the time response of the system under study. 

The availability of the intersar^le response is also an advan- 
tage to the designer, so that a design may be accomplished which 
accounts for the total time response and not that solely of the sam- 
pling instants . Therefore , the time domin matrix method has strong 
advantages over the commonly used z- transform on both of these points. 

Since adaptive systems may be treated as time varying systems, 
the time domain matrix method may be applied in general. Further- 
more, a learning control system with a predictor may be realized 
and investigated ty the use of tii® domain matrices , The investiga- 
tion of conditional control systems incorporating a model was also 
accomplished . 

For nonlinear systems, the time domain matrices provide a useful 
analysis and design technique. Ute characteristic of a deliberately 
introduced nonlinear compensator may be determined and investigated^ 
that is, it may be specified by the design approach. Furthermore, 
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the derivatives of the error or any other desired signal may be 
evaluated readily by matrix methods and used to determine the phase 
space response to aid the designer in the investigation. The use of 
a fictitous sampler and hold in a ccn^onuous data system permits 
these techniques of investigatiori and design to be applied. 

The design of the discrete coJET'ensator by means of the time 
domain matrix method provides the designer direct control of the 
time response of the control system. Therefore, the designer may 
specify the required system response and determine directly, the 
necessary discrete compensator. The i e.ilizability, stability, ar*d 
sensitivity of the discrete comrensatcr v>a3 Investigated, to provide 
the designer with an insight into the general characteristics of the 
discrete compensator. 

In comparison with other mothods cf investigation and design, 
the time domain matrix method is more <?,ccurate due to its numerical 
formulation. Furthermore, the accurac5>’ arid ease of calculation for 
the matrix method remains the same regardless of the order of the 
system transfer function, while the accuracy decreases and the dif- 
ficulty of calculation increases for other methods, such as the z- 
transform. The time necessary for calculation of the time response 
of a control system is considerably less than that of standard com- 
plex frequency techniques. Actually o the time response nay be de- 
termined as rapidly as the determination of instability or stability 
by standard techniques such as Routh's method. 

Tlie transfer matrix of the system may be easily determined by 
experimental pulse techniques tc sr. accuracy of at least three 
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percent. This method of system characterization may be used for the 
identification process in adaptive systems. 

Therefore, it has been shown that the time domain matrix method 
is applicable to a wide range of control system problems. Techniques 
have been developed which permit the solution of problems which are 
not solveable by any other method , Examples of such problems are the 
analysis and design of systems containing more than one nonlinear 
element, and the analysis and design of time varying systems, for con- 
tinuous or sampled data, 

A method has been developed which permits the design of a com- 
pensator for continuous or sampled data systems, without the use of 
trial and error methods. It also has been determined, that the 
difficulty of the application of this method and the labor involved, 
is considerably less than for other known methods when applied to 
systems of reasonable complexity. The advantages inherent in this 
method are sufficiently great, that it should find wide application 
in the engineering analysis and design of systems, 

?-3. Future Work 

The time domain matrix method may be applied to the complete 
spectrum of control problems. Therefore, the possibilities for future 
investigation utilizing this method are boundless , 

For all classes of control systems, it would be valuable to in- 
vestigate the problem of the location of the discrete compensator in 
the control loops, and determine and classify the relocation of the 
compensator on the time performance indices. 

For nonlinear systems, the investigation of the use of the error 
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derivatives and the signals available at intermediate points through- 
out the loops, would be very worthwhile. 

For time varying and adaptive system, this method would bring a 
new insight into the very difficult problems. It is possible to use 
the time domain matrix approach to investigate the method of steepest 
descent for adaptive systems, as one example. 
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APrEMDIZ A 



THE SYSTEM TRANSFER MATRIX 

A-1. Determination of the System Transfer I'ktrix Using the Z-Transform, 
As discussed in chapter two,, the determination of the system 
transfer matrix G is an important first step in the use of the time 
domain matrix method. It is necessary to evaluate the values of 
the values of the impulse response at the sampling instants. If the 
system transfer function is given in the laiplace variable, it is pos- 
sible, by standard methods, to transform to the z complex variable 
sT 

where z = e , Then the values of the impulse response are obtained 
by inverting the z-cquation to the time domain. The simplest method 
to accomplish this is use of division to yield the response at the 
sampling instants. That is, G(z) may be expanded as ° 

G(z) = go + gi7. (a-l) 

This series is the constant term and princip<al part of the 
Laurent expansion of G(z), , 

As an example , consider a system as shown in Figure A-l , 




Figure A-l. Open Loop Sampled System 
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The z transform of this system may be written as ° 



0h0(.) = ^ 



Z-1 

Then dividing we have 



s*(s l) ® 

K(„?679z“^ > o2642z“^) 



(a-2) 



U°e J 

z=o3679 



1 “ 1.3679z"^ + .3679z"^ 



=>3 

— = .3679z ^ •«“ ,7675s ^ ,9145z ^ + ••• 



(a-3) 



This method is tedious in that one must transform into the z variable 

and then divide. It also involves an inaceOTacy introduced by the 

2 

round-off in division. Jury presents a method of inversion avoiding 
division which uses G(z) in the form" 



G(z) = 



P2« 



o o C P ^ 



-n 



Po + Pi2 



- =2 =2 --m 

1 + q^z +9^2 •«» ... q^,,? 



(a -4) 



where m ^ n. 

Then , it can be shown that ° 



go 


= Po 






g1 


= Pi 


“ 9igo 


(a-5) 


g2 


= P2 


" qi gi “ 92go 






= Pn 


“ 9i gn- ■« " 9.2gn=-2 “ • ‘ ‘ 





Jury's method reduces the inaccuracies introduced by the division 
process , but remains inaccurate for a value of n greater than four 
due to the roundoff error of multiplication. For a higher order sys- 
tem, the use of the z-transform and inversion is unwieldly and inac- 



curate 



• • ^ 








I 



I 



I 










I 
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A-2. The use of the Impulse Response to Evaluate the System Transfer 
’Matrix . 

A simple and more direct method of evaluation of the S 3 rstem 
transfer matrix is the use of the impulse response of the system. 

If the impulse response g(t) is determined, then the elements of 
the system matrix are found by substituting t = nT, That is, 



g„ = g(t)| = g(nT) (a-6) 

I t=nT 

Since, invariably the impulse response is made up of a sum of 

time functions, the arithmetic operation involved in the evaluation 

is addition instead of division or multiplication. Therefore, the 

roundoff error of calculation is reduced as well as the difficulty 

of manipulation. Now reconsider the exam.ple of the previous section. 

The system transfer function is ° 

“ST 

G.G(s) = - (a_7) 

s^Cs+l) 

and the impulse response is I 



g(t) 
-sT 









s (s+1) ) 

The e simply implies a delas’’ cf cne sampling period. Then 
the inverse ist 

-ts 



g(t) = 



t -(1 



) J » j^(t - T) - (1 - e"^^ " 



] u(t-l 



-T) 

where T = 1 second as on figure A-1. (a-9) 

Then the values of interest at the sampling instants when t = nT = n, 
are* 



g = T + e""^ - 
“n 



1 + e 



e“(n “ 1 



(a-10) 
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Then, it is found’ 



go = 0 

gi = = o 36788 (a-11) 

g 2 = 1 + e”^ - = ,76746, etc. 

Therefore, the values of the sjrstem matrix are rapidly and 
accurately evaluated using a mathematics table . This calculation is 
a great improvement in accuracy over that of division in equation a-3. 

A-3. Evaluation of the System Transfer Jfetrix for Intersample Response 
Using the fbdified-z^Transformo 

In order to determine the response of a system between the sam- 
pling instants, it is necessary to obtain the system transfer matrix 
for values of time between the samples. In order to accomplish this 
with standard notation, let ra = 1 “ h, where A is the percentage de- 
lay from the sampling instants as shown in figure A-2. 




Figure A-2. Delayed Sampling 
Then, writing the transformation t (n - a)T one hast 
t = (n-A)T= (n-1 + n)T 0< m< 1 

and the modified z transform 1st 




i 



G(z,m) = 2 



(a-12) 



cc- 




k=0 



g kT •«= mT)z 



=k 



The output response of a system between the samples Y(z,m) may then 
be determined for an imput X(z) as I 

Y(z,ra) = G(z,ffi) X(z) (a-13) 

G(z,m) may be determined and inverted by division or an alternate 
method to obtain the elements of the system matrix for intersampling 
response . Then the intersample output response is I 

Y(m)] = [G(m)] xj (a-U) 



where 




go(m) 0 0 

g, (m) gc(m) 0 

g2 (ra) Si (ra) go (e) “ 0 » 



Since G(z) = zG(z,m) (a-15) 

m — ► 0 

the G ] matrix may be determined by allowing m-* 0. 

As an example of the determination of the modified system matrix 
using the modified z-transf orm consider the previous example . The 
modified z transform of the system is figure A-2 is” 



= ( m^-e“^^l )z^ ^ P.Z ^ (,3679m^ e~”^ - .7358) 

z^ - 1.36792 o3679 

where pi = 2,3679 - ra (1.36 79) “2c ™ 

Expansion by the method of section A=1 yields 

go (m) = Po (m) = m + e“‘“ - 2 

gi (m) = p, (m) - po (m)q^ 

gjCm) = • pjCm) - p, (m)q^ = Pe (°) (g? “ Gi ') 

and so on. 



(a-16) 



(a-17) 
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A-4. Evaluation of the System iTansfer ^&,trlx for Intersaunple Res- 
ponse Using the Impulse Response, 

The values of the modified mtrix are simple to evaluate by the 
impulse response method. One has“ 



g (m) = g(nT + AT) = g{t) j 

1 1 = nT AT 

= g(nT + T - mT) 

Then, for the example previously discussed one has I 



(a~lS) 



go(m) = (t - 1 + e ) 
= AT - 1 + e“^^ 



t - riT AT , T = 1 



(a-19) 



= A - 1 + c 



A = 1 - m 



Then, for example when A = 1/4 9 m - 3/4 one has° 
go (3/4)= - 3/4 = ,0288 



For n > 1 one has 



= T * - e" 



, T = 1 



t = nT + AT 



Therefore , 

g, vm) = 1 + 

g,w = 1 . 

and for m = 1/2 , A “ 1/2 

g^ (1/2) = .61660 = El 
g2(l/2) = .85R95 = g. ^ 

^ o ^ 



(a-20) 



(a-21) 



(a-22) 



where e . is another form of notation, 

®n+A 

A-5. The System I'tatrix for A Type I Second Order System with A Hold, 
The system transfer function for a type I second order system 
with a hold may be written as “ 

0 G(S) = (a-23) 

^ s^(s + a) 

Then using the inverse Laplace transform the impulse response is ' 



g(t) = ^at - ( 1 - for t ^ T 

a V 



(a-24 ) 



and 



P(t) = 



2 






)) - (a( 



at - = (a(t-T) - (l-e'^^^"'^) 



=at =a(t»T) 

aT + e “ e 



^ t>T 
(a-25) 



Therefore 



^n= T 

a 



j^aT + 



=anT =‘3T(n-l 

e “ e 



n 1 



(a-26) 



and 






A 

2 



a 



[aT 



=aT (n-s-A ) =aT (n+A~l ) 
e “• a 



(a-27) 



The matrix values for this second order system are tabulated in Table 
A-1 for three values of the constant aT. These matrix values my be 
used for any second order sj-^stem with a hold which possesses the 
sane aT value, making the proper allowance for the 1/^ factor ap- 
pearing in equation a-26. For example, a second order system with a 
pole at 10, a sampling period of .10 second, and a gain of 100 will 
possess the mtrix given in Table A=1 for the aT = 1 system with a 
gain of one . 
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A-6. The System I^trix for a II Second Order System With A 

Hold Network. 

The transfer function of a type two system of interest nay be 
written as " 



O^G(s) = 

h 3 

s 



Then using the inverse Laplace transform one obtains ; 

v2 



(a-?8) 



G(t) = K I (t + |“t^) 
Then, at the sampling instants. 



(t=T) l^t-T)' 



J u(t-T)| 



g = K i T + ~ 
®n 1 2 



[(nT)^ » ^(n-l)T^^]| n^ 1 



(a-29) 



(a-30) 



The values for the matrix of this system are given in Table A-1 , 
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TABLE A=1 



ELEIENT VALUES FOR THE SYSTEM TRANSFER ’^ATRIX FOR V.\RIOUS SYSTEl'^ 



la Typ« I Second Order System With A Hold Network 
GhG(s) ^ 

^ r = 1 second , aT = 1 



n 




n->A 


®(n+A) 


0 


0 


.5 


.10653 


1 


,3^7BB 


1.5 


.61660 


2 


.76746 


2.5 


.85895 


3 


.91445 


3.5 


.94812 


4 


.96853 


4.5 


.98091 


5 


.98842 


5.5 


.99298 


6 


.99574 


6.5 


.09741 


7 

? 


.99843 

.99943 

1.00000 


7.5 


.99905 



(n+A) 


^n+A 


(n-o-A) 


^n+A 


.25 


.0288 


.75 


.22237 


1.25 


.5077 


1.75 


.7014 


2.25 


.8189 


2.75 


,89016 


3.25 


.9334 


3.75 


.9596 


4.25 


.975 


4,75 


.9852 



Type I Second Order System With A 
GhG(s) 


Hold Network 


K 


s^(s + l) 


T = 


= 2 seconds , 


n 


®n 




^n>A 


0 


0 


c5 


.36788 


1 


1.13534 


1.5 


1.6819 


2 


1.8830 


2.5 


1.95695 


3 


1 .98416 


3.5 


1.99417 


4 

5 


1 .99786 
1.99971 
2,0000 


4.5 


1.99921 
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Type I Second Order System with a Hold Network 



-sT 

= — 5 T = .5 seconds , aT 

^ s2(s^-l) 



n 


^n 






0 


0 


,5 


.0288 


1 


.10653 


1.5 


.19357 


2 


.26137 


2o5 


.31413 


3 


.35523 


3.5 


.38727 


4 


.41221 


4.5 


.43163 


5 


.44674 


5.5 


.45853 


6 


.46771 






7 


.48041 






8 


.48812 






9 


.49279 






10 


.49563 







II A Type I Third Order System With A Hold Network 



\g(s) 

K 



s^(s+l)(s+2) 



T = 1 second 



n 


Sn 


n-^A 


^n+A 


0 


0 


.5 


.05824 


1 


.33616 


lo5 


.7845 


2 


1.18686 


2.5 


1.47885 


3 


1.67364 


3.5 


1.7983 


4 


1.87626 






5 


1.9540 






6 


1.9830 






7 


1.99372 






OO 


2. 0000 
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Type II Second Order ^stem V/ith A Hold Network 



Gj^G(s) 

if- 






T = .05 seconds 



g 

n 

0 

.9375 

1.3125 

1.6875 

2.0625 

2.4375 

2.8125 



Type II Second Order System With A Hold Network 



G^G(s) 

if" 



15(s4-10)g-e° 

_3 



11 



T = .10 seconds 



^n 


n-o-ii 


^n-t-A 


0 


. ^ 


.9375 


2.25 


1.5 


3.000 


3.75 


2.5 


4.500 


5.25 


3.5 


6.000 


6.75 






8.25 







Vl ^1*50 
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Figure A>3. The Impulse Response For Three Systems . 
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'iFFEMJn 3 



A IKTHOD OF II'nCRTIIC- A LC'^-R TFiTANGUL.iR I-IiTRIX 



The inversion of a lower triangular Jatrix is a necessary step 
in the evaluation of the response of a system by means of time do- 
main matrices. It was shown in section 2-6 that the output time 
response was given by equation 

(b-1) 



cj =[[!]- [[i]* [a]] -^| r] 



Therefore, it is usually necessary to calculate the inverse of 
[l] <• [o] =[a]. If [.ij is VT<*.t5n: 



[A] = 



1 



*^21 



>'22 



0 

0 



^31 ^32 ^33 



Then a mtrix equation nay be wrltten°'' 

[a] [&]-i y] = h] 

The solution for Y is° 

v] = [b] h] 

where [b] = | [a] [ b] | 

By the rules of matrix iuversiono 

[B]=[b][Aj-l 



or 



[a]-i =[6]-Mb] 



(b-2) 



(b-3) 



(b-4) 

(b-5) 



The equation b-2 may written fcr clarity as! 







• •* *■«>!*» 



I » r 

• ill • • 

•*'iif Ji)i . ( • <■( 
• j«jr !• 



Iff I ■*• 



7i 




( ’^)V^ 

3-1 1 


^ 72 


a 31 

(^)7i 

3-11 


a 

+ 

a^ 


etc , 





= h. 



(l>-6 ) 



3? 



22 



Then it can be seen that 
3-11 0 



72 



73 



h, 



[b] = 



^^22 



0 

0 



0 

0 

^33 

0 



0 

0 « 
0 

^44 



a 



nn 



and 



[6r^ = 



(l/a. . ) 
0 
0 



0 



(l/a,.) C 



0 



0 



11/ a 



as-* 



0 0 

0 0 

Q o o 






Hence, |]a] can be obtained by Biultiplloation of the successive 
rows of [b] by (l/a,,), (l/aj-i)^ etc, Tnen it can be seen that 

the elements in the first cclum:* of [bJ are the values of when 



H] = 



1 

0 

0 



= f 



5 O9 Op « ° « 



f 



The elements in the second column of j^B] are the values of y] when 
= " 1 ^ 0 , 1 , 0 , 0 , « ® *1 and so on,. 



1^*5 




4 



♦ 



I- 



I 




Ur 









Vi^ 
















% 




9 4 * 



• r * |»' -iMoi'lim 




When 



Hj = {l, 0, 0, 
7i = 1 



•1 



one obtains the first column as' 




(b-7) 



This solution suggests a method of computation which allows these 
calculations to be easily carried outo The elements of hy 

which the rows of j^Bj are to be multiplied, are recorded on the ex- 
treme right , This useful scheme , in the form of a table , is shown 
below for a 4x4 matrix. 




a 



31 



- (-^) - 



a 



a 



’’22 



10 0 

b,, 1 0 

^31 ^32 ^ 



0 

0 



1_ 

i 



1 - 

^22 



0 



1 



^33 



h 




- 




^•22 


3,33 






1 



a 



44 



The rules for setting up this table are. 

i) To form the left-hand array enter blanks in, and to the 
right of the principal diagonal. Derive the remaining 
elements from the A elements as shown, 
ii) Commence the right hand array by entering units in the 
principal diagonal and ze^os to the right of the diagonal 



iii) Calculate the remaining elements of [b] in succession 



by the following roethcd . To obtain b^^ postmultiply 
the row of the left-hand array level with b^j by the 
part colur,in of [bJ standing above j • Blank elements 
of the left-hand array are to be disregarded, 
iv) To find multiply the rows of [b] by the scalar 

factors on the right. 

For example, the first col'oian of [b] is completed below' 



■'21 



^31 



41 





^21 

. ” ^11- 


hi 








3' 31 






• —1 

1 




aTi” ’ 


^22 . 


] 


1^21 


= 


ail ’ 


„ ^42 

Q 

it 2 J* 


V , 





(b-9) 



^21 



^31 



As an illustration of this method consider the inversion of the ma- 
trix ^ A J given as I 



[a] = 



1 

0 

3 

1 



0 

-1 

-11 

-6 



0 

0 

-36 

-14 



0 

0 

0 

22 



(b-10) 



One then sets up the computation table following the given rules, 



— 


— 








0 


0 


0 


1 


0 


— 








1 


0 


0 


-1 


-3 


-11 






j *^31 


^32 


1 


0 


1 

“ 36 


-1 


-6 


7 

“ 18 




1 

H 


ha 


h3 


1 


9 

22 
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Then, the 



are calculated as fcllcvs 




bji 



^31 



= [0] [i] = 0 





^ 1 


[-3 , -11] 


'^211 







1 




= -3 




b 



32 





=11 



r 



0 



1 








77 

18 



18 




7_ 

18 



Then, one obtains 



[B] 




0 




0 

0 

1 

7 

18 



0 

C' 

0 

n ! 

JL. ! 



and since 




= [b]-MB] 



ori'& h,as° 




1 0 

0 -1 

i. H 

12 36 



0 

0 

36 



0 

0 

0 



(b-11) 



(b«l2) 



(b-13) 



For a physical system it was shewn ir* appendix A that invariably 



the element gg of the system matrix has a value of zero. For these 
physical systems one is usually interested in evaluating the inverse 
of [l] + [gJ and therefore the matrix i6j~^ is simply equal to[l]o 
As an illustration consider a type I system with a hold that has 



the following transfer function 

=sTh 



GvG(s) = p 



where K = 1 and T = 1 



It was shown in appendix A that the system matrix ist 



[a]= 



i 0 


0 


0 


0 


0 


0 


.3679 


0 


0 


0 


0 


0 


.7675 


.3679 


0 


0 


0 


0 


.9145 


.7675 


.3679 


0 


0 


0 


to 

• 


.9145 


.767? 


.3679 


r. 


0 


.9884 


.9685 


o 9 U 5 


.76"’5 


.3679 


0 



(b-U) 



Then, one obtains for ] •«• [gJ j, the principal diagonal containing 
element values of unity. Completing the table for inversion, one has! 



— 
















-.3679 


_» 














-.7675 


-.3679 


— 












-.9145 


-.7675 


-.3679 












-.9685 


-.9145 


-.7675 


=.3679 










-.9884 


-.9685 


-.9145 


-.7675 


-.>679 










1 


0 


0 


0 


0 


0 


1 




-.3679 


1 


0 


0 


0 


0 


1 




-.6321 


-.3679 


1 


0 


0 


0 


1 




-.400 


-.6321 


-.3679 




0 


0 


1 




+ .001 


o 

o 

• 

1 


-.6341 


-.3679 


1 


0 


1 




+.2528 


+ .001 


-.-.oc 


-.6321 


-.3679 


1 


1 



1 c 



Furthermore , for non-time v-ir/ing systens the [g] mtrix always 
contains equal values along any uprf;: ,i:ft to lower right diagonal. 
Therefore, inversion of [l] ^ [u] need be accomplished for only one 
column as is obvious from the example , This class of mtrix called 
a diagonally invariant mtrix (D„I„Mc) in chapter three, lends it- 
self to the use of a recursion form'ola for the elements of the in- 
verse matrix. For the first ccl*ur':n of the matrix [b] one may write* 



(s 



a. 



" ^21 ’=k2 ’ 



(b-15) 



"kl '"kl "k-l,.l i 1 

For example, in the previous solution* 

bj, — — (a^, bjj + ^21^32) (b- 16 ) 

Since, the elements are invariant on the diagonal, one has* 

^3 3 ~ 1 ~ ^ 



(b-17) 



and 



^32 ~ ^21 



Therefore , 



b3i = -(aji + ajib 



21 / 



U7675 „ 3679 (-.3679) ) 



06321 



In order to obtain accurate resu3.ts, it is necessary to use 
either a desk calculating machine or an automatic digital computer. 
This is necessary in order to avoid the roundoff errors of multiplica- 
tion, For the usual requirements of one percent accuracy, the desk 
calculator is completely adequate . 
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